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ABSTRACT 


Aloioat whole of the year, mid- to hifih-lattitude 
regions are covered with snow. The snowpack undergoes several 
cycles of melting and refreezing throughout the year, which 
influences the weather conditions, and surface and subsurface 
water. The capability of microwaves to penetrate snow layer and 
respond to variations in subsurface properties along with 
all-weather and all-time observational capabilities make Microwave 
ReiJiote Sensing an useful tool to monitor snowpack properties on 
global scale. The microwave portion of electromagnetic spectrum is 
advantageous additionally because of the large difference in the 
dielectric constant of liquid and frozen water, which causes a 

strong variation in the observed microwave signal when liquid water 
is present. 

In order to study the microwave remote sensing 
response over snow covered areas, detailed numerical calculations 
have been carried out over various two-layered and multi-layered 
models representative of such areas. The results show that the 

microwave remote sensing can be used to map the snow cover extent, 
to classify snow types, to estimate snow thickness, to detect onset 
of snowmelt, to estimate the free liquid water content of snow, 

to monitor the roughness of snow surfaces, and to estimate the lake 
ice thickness. The detailed discussion of the results have been 
presented in view of the proper selection of operating parameters 
of microwave remote sensors for the mapping of snow covered 

regions. 


(iii ) 



CHAPTER 1 


INTRODUCTION 


l.i CHSNERAL VIEW 

Th« volume of ice and snow on earth is enormous: 

enough to cover its entire surface with a layer roughly 60 
meter thick. Roughly 50 percent of earth’s surface is, at present, 
covered by snow or ice during most of the year. Ice is not 

uniformly distributed over the earth but is instead concentrated 
in the polar regions and mountains, attaining thicknesses of over 
4,000 meters; snow is concentrated in mid - to high - lattitude 
regions (EOS, V.II f). Figure 1.1 shows the sea ice, glaciers and 
snow regimes for the northern and southern hemisphere. 

Snow after being deposited on earth’s surface serves 
as a source of water vapour input to the atmosphere through the 
processes of sublimation and evaporation and also as a source 

of water to the soil and river systems after its melting. Uith 
the increasing demand for water throughout the world, effective 

management of the water resources is required in order to predict 
and distribute the available water supply efficiently. 

The use of remote sensing techniques offer an 
outstanding opportunity to survey the snow cover much more 
extensively than can be done practically from the ground, since 
snow fields normally occur in remote and inaccessible areas where 
the data collection becomes very difficult and hazardous. 




NORTHERN HEMISPHERE 



WINTER SUMMER 


SOUTHERN HEMISPHERE 


KEY 

sea ice 

jnPMI glaciers 

SNOW 


Fig. 1.1. Sea ice, glaciers and snow regimes for the northern and 
southern hemisphere (EOS, v.IIf). 
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1.2 CCFINITIOII AMD CCMCEPT OF REMOTE SEMSIM6 


Renottt sensing la the technique of obtaining 
inforaation about an object, area, or phenonenon through the 
analyaia of data acquired by a device that ia not in contact with 
the object, area, or phenoaenon under inveatigation. Reaote 
aenaing ia poaaible in three regiona of electroaagnetic apectrua 
(Figure 1.2) naaely, vlslbl* and n*ar Infra-rad, -tharmal infra-rad 
and alcrowava. The two baaic proceaaea involved are 'data 
acqulBltion' and ’data analysis’. The eleaenta of the data 
acquiaition proceaa are : 

(a) energy aourcea, 

(b) propagation of energy through the ataoaphere, 

(c) energy interactiona with earth aurface featurea, 

(d) re-tranaaiaaion of energy through the ataoaphere, 

(e) air borne and/or apaceborne aenaora, and 

(f) generation of aenaor data in pictorial and digital fora. 

The data analyaia proceaa involvea exaaining the 
pictorial and digital aenaor data. In short, we use aenaora 
to record variations in the way earth aurface features reflect 
and eait electroaagnetic energy, which is used further to deduce 
the aurface and subsurface inforaation. 

1.3 WHY MICROWAVES FOSt REMOTE SENSING 9 

The ndcrowavef portion of electroaagnetic apectrua 
extends froa 0.3 to 300 GHa range of frequency (1 aa to la range of 
wavelength). The aicrowave region is further classified in various 
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bands of narrow fraquancy range and ia shown in Tablet. 1 (Ulaby at 
al., 1981). Reaote sensing by alcrowaves has additional advantages 
over the visible and near infra-red and thernal infra-red 
regions of electroaiagnetic spectruai because - 

(a) nicrowaves have capability to penetrate clouds and to sone 
extent rain, i.e., all-weather observation is possible, 

(b) Microwave sensors have capability to probe the target day and 
night , and 

(c) Microwaves have capability to penetrate deeper into the 
target and consequently infornation with depth is available, 
whereas in visible and infra-red region only surface 
inforMation is known. 

1.4 CejECTIVES OF THE STUDY 

In developing countries like ours. Monitoring of snow 
and ice cover is very iMportant in water resources developaent and 
ManagSMent. By knowing the physical paraMeters of snowpack, an 
accurate estiMate of annual runoff of Major rivers due to Melting 
of glaciers (alMOSt all the rivers in northern part of India) can 
be predicted, that will aid the effective nanageaent of water 
resources to distribute the available water supply in donestic, 
irrigations! and industrial sectors. 

The present study was carried out with the following 

objectives : 

(i) to relate the dielectric constant of snow with frequency, 
density of snow, liquid water content and physical 
teaperature of the snowpack. 
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TABLE 1.1. Band daslgrutllons In the microwave region 


Band 

Frequency Range 

C&izy 

P 

0.225 

- 0.39 

L 

0.390 

- 1.55 

S 

1.550 

- 4.20 

C 

4.200 

- 5.75 

X 

5.750 

- 10.9 

K 

10.90 

- 36.0 

0 

36.00 

- 46.0 

V 

46.00 

- 56.0 

u 

56.00 

- 100.0 
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(il) 


to oatioiate the scattering, absorption, and extinction 
losses of electroaagnetic energy and penetration depth of 
aicrovave signal in the snowpack, 

(iii) to discrininate snow covered areas froa bare ground 
surfaces, 

(iv) to classify snow types, 

(v) to estiaate snow thickness for the deteraination of snow 
water equivalent, 

(vi) to estiaate the liquid water content of snow, 

(vii) to study the effect of roughness of snow surfaces caused due 
to snow aelting and local winds on aicrowave reaote sensing 
response, 

(viii) to study the behaviour of aicrowave signals over natural 
lake or river ice (in case of pure and saline water), 

(ix) to study the effect of incidence angle, polarisation (both 
horizontal and vertical) and frequency on reflectivity, 
eaissivity and brightness teaperature of snow and ice 
covered surfaces, and 

(x) to suggest ultiaately, the optiaua sensor paraaeters for 
reaote sensing of snow and ice cover. 

l.S ORGANIZATION OF THE STUDY 

Besides this introductory chapter, the thesis contains 
four aore chapters, covering the whole aspects of aicrowave 
reaote sensing to aonitor snow and ice cover on the surface of 
the earth. Chapter II describes the general theory of passive 
aicrowave reaote sensing systeas. A detailed discussion of 
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«ala«ion behaviour of two-layered and nultl-layered snow and ice 
surfaces has been dealt here. The effect of surface roushness 
on reflectivity caused due to local features like strong 

katabatic winds and snow nelting is also discussed here. Chapter 
III deals with the physical properties of snow and ice as well as 
various algorithas to calculate the dielectric constant of snow, 
ice and water. The behaviour of real and imaginary parts of 
dielectric constants of dry and wet snow with the frequency of 
microwave signal, snow density, liquid water content of snow and 
surface temperature has also been discussed in this chapter. In 
chapter IV, we have presented the microwave remote sensing response 
over various layered models representative of snow covered 
terrains. Conclusions and further recommendations are given in 
Chapter V. 

The Computer programs developed and used in the 
present study are appended at the end. 



CHAPTER 2 , 


THEORY OF MICROWAVE REMOTE SENSING 


a.l INTR<X>UCTION 

All-w«4ith«r and all-tl»« obaarvational capabilities 
and deeper penetration into the ground aake aicrowave eenaore 
nore advantageous over visible and infrared sensors. In aicrowave 
reaote sensing the aicrowave response (scattered or eaitted energy) 
is controlled by the largat as well as s*naor paraaeters. Target 
paraaeters include- i) surface geoaetry, ii) subsurface geoaetry, 
and iii) dielectric properties of the surface and subsurface 
aaterial; and sensor paraaeters include i) wavelength or frequency 
of electroaagnetic radiation, ii) polarization, iii) angle of 
observation, and iv) spatial resolution. The fundaaentals of 
aicrowave reaote sensing systeas as an overview and various 
foraulae to coapute the brightness teaperature over two-layered and 
aulti-layered aodels are given in this chapter. 

2.2 MICROWAVE REMOTE SENSING SYSTEMS 

The aicrowave reaote sensing systeas are classified 
into two categories viz. Passive and Active, based on the source 
of electroaagnetic energy. The systeas which detect the 
eaitted aicrowave radiation by the target are called passive 
systeas, e.g., aicrowave radioaeters. The systeas which transait 
the aicrowave radiation to illuainate the target and at the saae 



time receive the return or back scattered siftnal are called 

active microwave remote sensing systems, e.g. , scatterometers , 
altimeters, side looking airborne radar (SLAR) and synthetic 
aperture radar (SAR) etc. 

The microwave remote sensing systems can also be 
classified as non-imaging and imaging sensors. The examples of 
non-imaging sensors are radiometers in non-scanning mode, 
scatteromet er , altimeter etc. and are very useful for understanding 
the physical phenomenon and modelling purposes. In the case 
of imaging sensors, the information is gathered and the results are 
depicted in two-dimensional pictorial form known as images. 
The radiometers in the scanning mode, side looking airborne 
radar (SLAR) and synthetic aperture radar (SAR) are some of the 
the examples of imaging sensors. These sensors produce the image 
of the object, which is useful for the detection and assessment of 
natural resources. 

2.3 PASSIVE MICROWAVE REMOTE SEl^MG 

A microwave radiometer like other passive 
electromagnetic sensors measures the emitted energy from the 
targets within the field of view of its antenna. This 

emitted energy is quantified in terms of brightness temperature 
(Tg) which is a well defined function of electrical property, 
i.e., dielectric constant of the medium under investigation and 
also of physical temperature of the object. 

2. 3. 1 Derivation at Brightness Temperature 

According to Planck’s quantum theory of black body 
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radiation 


all tha objecta abov* O^K radiate electronaftnetic 
energy. The aatheoiatical foratula given by Planck is 

2h£^ 

Bf “ , (2.1) 

C (h£/kT) - 1 ] 

where 

—2 —1 —1 

ie the blackbody spectral brightness, U a Sr Hz , 
h is the Planck’s constant, d.63 * 10 J-sec, 
f is the frequency in Hz, 

-23 

k is Boltzaann constant, 1.38 x 10 J/K, 

T is absolute surface teaperature in *^K, 

A 

c is the velocity of light, 3 x 10 a/sec. 


In microwave region, hf/KT << 1 (Rayleigh - Jean’s 
law). Therefore, the Planck's formula reduces to 


B = k T . (2.2) 


All the materials, referred to as grey bodies emit 
less than a black body does and do not necessarily absorb all 
energy incident upon it, so 


b’ » — k T- , (2.3) 

c 

* 

where is the spectral brightness of grey body. The emissivity 

t 

ia defined as the ratio of the spectral brightness of the grey 
body to that of the black body. Therefore, 


i 



(2.4) 
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or Tg . * . T , (2.5) 

where e is referred to as enissivity of the nediun and varies 

between zero and one (0 < e < 1). Eaissivity of a Material 
depends upon frequency, polarization, incidence an£le, conplex 
dielectric constant, surface teaperature and surface roughness. 

2.3.2 BrlghlzMSs TAnqpMeratura Retcordsd by RadlomstArs 

In general, the radiation received by an airborne 
or spaceborne radiometer can be split into three components: the 

first is the self-emission from the surface of the object, the 

second is the reflection at the surface of the downwelling 
radiation from the atmosphere, and the third is the upwelling 
radiation of the atmosphere directly incident upon the radiometer 
(Figure 2.1). If the intervening atmosphere is lossy, the 
radiations emanating from the object will be attenuated and the 
lossy atmosphere itself will emit thermal radiation. Thus 
brightness temperature of the object can be written as 

X « T(rT, +(l-r)T ^)+T^ , (2.6) 

where r is the surface reflectivity and t is the atmospheric 
transmissivity. The first term in the . equation (2.6) is the 
reflected sky brightness temperature which depends on the 
wavelength and atmospheric conditions, the second term is the 
emission from the -surface of the object and the third is the 
contribution from the atmosphere between the surface and the 
receiver. 
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£Sm 3a 3 RA<li.oiMi''Ler’S 


A radioii«t«r ie ni«rely a paaaiva receiver which 
records self-emitted and/or reflected radiation by the scene. 
Host of the microwave radiometers are Dicke type of radiometers. 
A block diagram of the Dicke radiometer is shown in Figure 2.2 
(Ulaby et al.,1981). It consists of basically high gain antenna, a 
Dicke switch, a noise source, predetection section, square-law 
detector, synchronous demodulator placed in between the square-law 
detector and the low-pass filter (integrator). The predetection 
section consists of the RF amplifier, mixer, and IF amplifier and 
is characterized by a power gain (G) and bandwidth (B). 
Dicke switch, connected at the receiver input (at a point as close 
to the antenna as possible) is used to modulate the receiver input 
signal. The radiometer output is in millivolts. Calibration of 
radiometer gives the temperature corresponding to radiometer output 
voltage. 

3. 3. 4 Kasolutlon and Llstitations 

Two of the most important parameters of the 

radiometer are the ground resolution aind measurement accuracy. 

The specified temperature measurement accuracy of presently 

available radiometer is 0.1 ^K. The ground resolution of microwave 

radiometer depends upon the height of the platform and on the 

beamwidth of the radiometer antenna. From a platform flying at 

an altitude H, the linear extent of resolution will be 

X 

r “ H . = 1.2 H , (2.7) 

® D 
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wh«r« 

is ths antenna beanvidth, 

X is the wavelength of the receiving radiation, and 
D is the diameter of the antenna aperture. 

Froa the above equation, it is clear that for a 
practical antenna size, the achievable ground resolution fron the 
satellite height is poor. 

2.4, SNCtW PARAMETERS AMD SENSCMt RESPONSES 

Table 2.1 shows the snow paraaeters affecting aeasured 
response by various sensors. As shown in table, water equivalent is 
the key property for gasuaa ray sensors, albedo (the ratio of 
reflected to incoaing solar radiation) is aost iaportant for 
visible and near infra~red sensors, eaissivity is aost aeaningful 
for the theraal infra-red sensors, and the dielectric constant is 
aost significant for aicrowave sensors. 

Microwaves have the capability to penetrate snow and 
respond to variations in subsurface properties. The aicrowave 

sensors have potential to monitor aany iaportant snowpack 
properties such as depth, water equivalent, presence of aelting, 
and the state of underlying ground because emitted, reflected and 
scattered electromagnetic radiations from the snowpack are 
sensitive to changes in snowpack properties (Foster et al . , 1984). 

2.S MCH^L PARAMETERS OP SNOW 

The absorption and scattering coefficients are the 
model parameters (also known as propagation parameters), used to 
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TABLE B. 1. SSnow par Aiikat.*rs ng s*nsor r«sponB*s 


GaiMMa Ray Sensors 

- snow depth 

- background radiation . 

Vlslbl*.^tl*ar— Infrarad CAlb*do> Sansors 


- crystal size 

- contaainants 

- snow depth (only for shallow snow, i 
a few centiaeters) 

- liquid water 

- surface roughness 

Tharnal InfTarad Sensors 


- teaperature 
crystal size 

- liquid water 

Mlorowav* Sensors 


- liquid water 

- crystal size 

- water equivalent depth 

- stratification 

- snow surface roughness 

- density 

- teaperature 


e. upto 


soil condition 



coaputft the emission fro* a snow layer. The sun of the absorption 

(k ) and scattering coefficients (k ) Is known as the extinction 
® 0 

coefficient (k ) , which Is related as (ulaby et al . , 1986) 

w 

k. - • C2.8) 

The albedo (a) of the snow nedlua Is defined as 

a » k^ / k^ . (2.9) 


Scattering of alcrowave signals In snow Is aalnly due to the 
presence of Ice particles, therefore, scattering coefficient is the 
sua of the scattering cross-sections of all the ice spheres 
contained In a unit voluae, expressed as 


k 


s 



j«l 


(rj,n) , 


( 2 . 10 ) 


where is the nuaber density of ice spheres. The absorption 
coefficient of snow consists of two contributions: 1) absorption 
by ice spheres absorption by the background 
aedlua Therefore 


where 


"ai 



+ k 


ab 


» 


I 9. . 

j»l 


( 2 . 11 ) 


( 2 . 12 ) 


and 

^ab = 2kQ (1-vp n^', 

« 2kp (1-v^) jia ( ^^^)| . 


(2.13) 


In the above equations, k^ Is the wave nuaber In free 
space, v^ Is the voluae fraction of Ice particles In the alxture. 
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^ 1 ^ i« the iaaginary part of coMplex index of refraction of the 
background sfiediuai, and *|ji« the conplex dielectric constant of 
the background aiediua. The wave nunber (ICq) and volume fraction of 
ice particles (Vj^) are given by 


kg = 2fr / Xg , (2.14) 

, (2.15) 

where p^ and p^ are the density of snow and ice respectively. For 

9 f 

dry snow, the background medium is air and nj^ »0.0, hence k^= ^si ‘ 

For wet snow, k^^^ can be estimated by equation (2.13). The Rayleigh 

expressions Q (r,,n) and Q (r,,n) for an ice particle of radius r, 

are given by (Ulaby et al., 1981) 


8 n 




1. ( -K ) , 


and 




128 « , , 

— r I’'! 

3 X' ^ 


(2.15) 


(2.17) 


where 


and 


- 1 

K » 

jc + 2 




/ 


(2.18) 


(2.19) 


In the above equations, and are the complex 
dielectric constant of ice and background medium respectively. From 
equations (2.10) and (2.12), we can write 


128 n' 


N 


I^I 


j ’ 


3 X. 


j = l 


( 2 . 20 ) 



( 2 . 21 ) 


k 


ai 


8 n' 


I« ( - K ) 



J»1 


r 


3 

j ‘ 


If th* aean or typical radius of ice particles is r 
and all the particles are assuaied to be of the sane size, then the 
nuttber density of ice particles is given by 


N 




( 2 . 22 ) 


In wet snow, background nedius consists of air and 
water. In the present calculations water droplets are assuaied to be 
spherical and air as the host dielectric naterial, then the 
dielectric constant of the background nediun can be represented 
by Polder'-Van Santen two-phase nixing foraula (Ulaby et al., 1986) 


= 1 + 


(^ - 1 ) 
' w ' 


I 


u»a,b,c 


1 + A (jf> Its. 
u w b 


1 ) 


(2.23) 


where a is the voluaetric water content 

V 

o£ 

snow. 


is 

the 

dielectric constant of water, and 



and 

^c 

are 

the 

depolarization factors of spherical 

water 

droplets 

(- 

1/3). 

On 


solving, equation (2.23) sinplifies to 


2 + tfj^((^^-l)(l-3*^) - 1) - = 0 , (2.24) 

which is a sinple quadratic equation that can be solved to 
deteraiine the dielectric constant of the background nediua. Since 
ice is a scatter-free nediun, therefore total loss of 
electroaagnetic energy is only due to absorption which is expressed 
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aa 


k 

a 


kg 


/ 


) 


(2.25) 


wh«r« and are the real and iaaginary parts of conplex 

dielectric constant of pure ice. The penetration depth (*^p) 
aicrovave signal in a scattering aediua (like snow) is given by 

» 1 / , (2. 2d) 

where k is the extinction coefficient of snow layer. 

fS 


2.6 APF>ROAC3i TO COMPUTE BRl WITNESS TEMPERATURE 

In passive microwave remote sensing, basically two 
approaches are possible: the first is the non~coh*r*nt approach in 
which the theory is developed considering the basis of 
ttlaclromagna-tic Intanslby, and second is the cobarant approach in 
which alacbromagnatlc fiald is considered. In the present work, 
brightness temperature has been computed using non-coherent 
approach given by Ulaby et al . (1981) for two- layered case (as 
shown in Figure 2.3). The brightness temperature (Tg) at a point 
above the boundary between air and a medium consisting of two 
layers can be determined by the following expression - 


1 - r 


1 - r. 




*^1 ‘‘2 


r [(■ • t1(’ ■ t-K' - •]'. • M'=l 


1 - 




(2.27) 


where r^^ and are the Fresnel’s reflectivities at the two 

interfaces, L and a are the loss factor and albedo of the second 
layer respectively, and Tj^ snd T 2 are the average temperatures of 
the layer 1 and 2 respectively. The loss factor of the first layer 



is fiiv«n by - 


L * mxp ( k d 8 «c ©- ) , (2.28) 
wh«r« ©2 1 k^ and d reprsssnt th* refraction ansle, extinction 
coefficient and thickness of the layer 1 respectively. If the 
average physical temperature of the layer 1 and 2 are assumed to be 
the same (Tj=T 2 *T), equation (2.26) reduces to 

Tg - e . T , (2.29) 
where e is the incoherent emissivity and can be represented as - 


(1 - r^) 



[[■ • t](' - t1(- - ■) • (-^J) 


(2.30) 


The normalized impedances of the three media are 
given by the following expression - 


{ 77 ^ cos ©^ , for vertical polarisation 
77 j^ sec ©^ , for horizontal polarization ’ 

where 

T?i “ 1/ Vi ! . 

and i s 1,2,3. 


(2.31) 


(2.32) 


Here, and ^ . denote the relative permeability and 
relative complex dielectric constant of the respective medium. The 
refraction angles ©2 and ©^ can be determined by using Snell's law. 
The reflectivity at the two interfaces is given by - 

rj - ICZ 2 ” ^1^ / ^^2 ^1^1^ ’ (2.33) 

^2 * I ^^3 " ^2^ ^ ^^3 ^2^1^ * 




(2.34) 



As dJlscusssd in the next chapter, anow la deposited in 
discrete layers and the characteristics of one layer differ froai 
another, i.e., snowpack parameters vary with depth. The dielectric 
constant also varies accordinfily, which affects the emission 
behaviour of the snowpack. The brightness temperature for 
multi-layered case can be computed in the similar manner. The 
formulation discussed above was used to study the behaviour of 
brightness temperature over two-layered and multi-layered snow 
models. 

2.7 EFFECT OF SURFACE ROUGHNESS 

The polar and mid- to high-lattitude regions 
are surrounded by low pressure and the surface topography of the 
snow and ice surfaces are often disturbed due to the presence of 
snow melting and strong katabatic winds. The microwave signals 
incident upon these disturbed surfaces are partly reflected in the 
specular direction and partly scattered in all directins. In order 
to study the effect of surface roughness on reflectivity, we have 
used the following semi-empirical relation (Choudhury et al . , 1979) 

9 

r ■ r*** e ** cos^ & , (2.35) 

where r**^ is the specular surface reflectivity, & is the angle of 

9 

incidence and h is surface roughness parameter given by - 

h’ - 4 kj , (2,36) 

where k^ is the wave number in free space and o* is the surface 
height variation. 
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CHAPTER 3 


PHYSICAL AND DIELECTRIC PROPERTIES OF SNOW AND ICE 


3. 1 IMTRCXMICTION 

Th« ttlectrofiA£netic propertlea of & snov/pack vary 
aa its structura and conposition changa. As discussad in tha 
pravious chaptar> nicrowava sansor rasponsa is highly dapandant upon 
physical and dialactric propartias of a snowpack. An accurata 
analysis of tha nicrowava ranota sansing data raquiras prior 
knowladga of physical and dialactric propartias of watar, snow, ica 
and othar gaologic natarials (Singh,. 1984). In tha following 
sactions physical and dialactric propartias of snow and ica, and 
anpirical ralations for calculating tha dialactric constant of 
watar (both pura and salina) and bar a ground are discussad. 

3. a PHYSICAL PRC»>ERTIES COF SNOW AND ICE 

In araas whara annual snowfall axcaads tha loss of 
snow dua to nalting, a layar of snow gats accunulatad every year. 
Tha yearly accunulation of snow acquire innensa size and shape. 
Tha newly fallen snow, which is extraaaly light and highly porous 
nass of frozen watar and air, undergoes a series of changes with 
tiaia. It is first conpactad to a granular mass of alnost spherical 
grains with specific gravity increasing from as low as 0.05 to 0.50 
and porosity reducing fron original 95 percent to 50 percent or 



1«S8. This ia the firn or ne've’, a variety of ice. With 

continued conpaction due to increasing mass of snow layers from 
above, firn undergoes further changes towards a very coapact , 
dense, coarsely crystalline aass with a specific gravity of 0.916 
and with little or no porosity, defined as ice. 

Snow is a fine grained material with a large 

specific area. The shape and size of the snow particles change 
very rapidly during melting and refreezing cycles, especially 
in response to changing weather conditions. The snowpack 

properties vary throughout the life history of a snow cover and 
since snow is deposited in discrete layers, the pack may vary 
throughout its depth. Moreover, the effects of wind and 

local melting can create horizontal inhomogeneity. This 

variability makes monitoring snowpacks all the more challenging 
and crucial (Foster et al . , 1984). 

Fresh, dry snowflakes display crystal facets because 
they grow rapidly in the atmosphere. The flakes metamorphose, 
either via destructive metamorphism to rounded grains or via 
constructive metamorphism to large, faceted crystals known as 
depth hoar. 

Equilibrium thermodynamics favours destructive 
metamorphism, which controls the shape of the small, well rounded 
grains that generally dominate a dry snow cover. But when there 
is a steep temperature gradient, constructive metamorphism occurs, 
and large, angular, and hollow crystals of depth hoar grow in 
the warmer layers at the base of the snowpack (Foster et al., 
1984) . 

When snow is wet, its structure and its 
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elttctronagn«tic behaviour change Markedly. Two geometries of wet 
snow are coanon: freely draining snow with low liquid water 
content by volune) and snow with a high liquid water 

content where water flow is impeded. The grain size increases 
rapidly following even slight wetting. At higher liquid water 
contents, the grains grow at the faster rate than in freely 
draining snow. This type of snow is cohesionless, with well 
rounded grains about 1 mm in radius, and fully surrounded by water. 
(Colbeck, 1982). 

3.3 EOUATIONS FOSt CALCULATING DIELECTRIC CC»^ANT OF ICE, SNOW, 
WATER AND BARE <»OUND 

In microwave remote sensing, knowledge of the 

dielectric properties of the materials is important for the 

interpretation of the microwave reflection or emission data (Singh, 

1984). The dielectric properties of the materials are characterized 

* 

by its relative complex dielectric constant ( ) and loss tangent 

(tan <5^), where the subscripts r and s refer to relative and snow 
respectively. However, for simplicity the use of the term relative 
and hence the subscript r is avoided in the forthcoming 

fc 

discussions. The dielectric constant (« ) is a measure of the 
response of the material to an applied electromagnetic field. This 
response can be split into two factors: the first determines the 
propagation characteristics of the wave in the material, i.e., 
velocity and wavelength, and the second is a measure of the 
energy losses in the media. These two factors are represented by 
the r«al ( 4 ? ) and imaginary (« ) parts of the complex dielectric 

constant. The ratio 4?”/ 4> is called the loam tangant (tan 6^} for 
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the aateriAl. 


* • » » 


19 f 

tan <5 ® jc f £ 

8 ' 


C3.1) 

(3.2) 


where i 
pttrnd.tt.lvl'ty 


/ 

of 


-1 , £ ia 

the jMediun and £ 


» » 


referred to aa the dlAlactrlc 
la known aa the dlAl«ctrlc loss 

9 9 9 


factor of the nediua. In general, the valuea of £ and £ are 
functiona of frequency and tenperature both. 

The equationa for calculating the dielectric 
conatant of pure ice, aaline ice, anow, pure water, aaline water 
and bare ground have been diacuaaed in the following aectiona 


3. 3. 1 Dleloctrlc Constant of Pure Ice 

The real part of conplex dielectric conatant of pure 

f 

ice ( £ ) is Independent of both teeyMrature and frequency in 
the microwave region, and ia assigned the conatant value 

= 3.15 , (3.3) 

f 9 

however, the imaginary part ( ) exhibits strong variations with 

9 9 

both parameters. Valuea of are taken for the calculation of 

dielectric constant of anow and saline ice as shown in Figure 3.1 
at -1 and -20®C (Stiles et al . , 1981). The subscripts i refers to 
pure ice. 


3. 3. 2 Dloloctrlc Constant of Salin* Ice 

Saline ice ia a much more complicated medium, both 
structurally and electromagnetically than fresh-water ice and is a 
heterogeneous mixture of liquid-brine inclusions and air pockets 
interspersed within the ice medium. The brine inclusions. 
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Loss Factor of Pure Ice 



0.1 0.3 1 3 10 30 100 

Frequency (GHz) 


Fig. 3.1. The imaginary part of complex dielectric constant (jc. ) 
of pure ice. 
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containing salt and watar, show strong influence on the conplex 
dielectric constant of the saterials (Rankin and Singh, 1985, and 
Singh et al . , 1989) because of their high conplex dielectric 

constant when compared to that of ice. 

In general, the complex dielectric constant of saline 
ice i« a function of the following parameters - 

(a) the complex dielectric constant of pure ice 

(b) the conplex dielectric constant of the brine pockets or 
inclusions 

(c) the fraction of brine by volume (Vj^), and 

(d) the shape and orientation of the brine pockets or inclusions. 

The salinity of liquid brine is governed by its 

temperature and can be calculated by the following empirical 
expressions (Assur, 1960 ; Poe et al., 1972) - 

S. » 1.725 - 18.756 T - 0.3964 , -8.2 :£ T i -20 (3.4) 

D 

S. « 57.041 - 9.929 T - 0.1620 - 0.002396 , 

D 

-22.9 ai T S -8.20 *C (3.5) 

a 242.94 + 1.5299 T + 0.0429 , -36.8 < T i -22.9® C 

(3.6) 

S,. » 508.18 + 14.535 T + 0.2018 , -43.2 ^ T S -36.8® C 

(3.7) 

The salinity of liquid brine is related to 

normality by (Klein and Swift, 1977) 

[1.707 X 10“^ + 1.205 x lO"^ S^^ + 4.058 x lO"’ S^ J, 

(3.8) 

» f f 

and and for pure NaCl solution are given by Stogryn’s 

D D 
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(1971) fornul&tion as 


bO wot 


4i. ■ * + 

b VKB 


1 + (2fi£Tj^)' 


(3.9) 




bO woo 


(1 + 2n£T^)‘ 


2n £jc. 


(3.10) 


where « » 4. 

woe 

9 , and 


^bO^^'**b> 

- (T,0) a^(N^^) , 

(3.11) 

% (T'V 

» (T,0) b^(T.Nj^) , 

(3.12) 

% (T.Nj^) 

= (25, Nj^) Cj^CA ,N^^) . 

(3.13) 

The above £unctions are given by the 

£ollowing 

expressions - 




88.045 - 0.4147 T + 6.295 x lO"^ 



+ 1.075 X lO"^ , 

(3.14) 

^ (Nfe) - 

1.0 - 0.255 Njj + 5.15 x 10“^ - 6.89 x 

-3 3 

10 , 



(3.15) 

^^^(1,0) » 

(1.1109 X 10“^° - 3.824 x lO"^^ T + 6.938 

X lO-^V 


-5.096 X lO"^*^ T^) / 2n , 

(3.16) 

bj(T,Nj^) - 

1.0 + 0.146 X lO"^ T - 4.89 xlO"^ 



- 2.97 X 10~^ + 5.64 x 10“^ , 

(3.17) 

%(25.Nj^) « 

Nj^(10.39 - 2.378 + 0.683 - 0.135 



+ 1.01 X lO"^ N* ) , 

(3.18) 

Cj(/1,N^^) » 

1.0 - 1.96 X lO"^ + 8.08 X 10~^A^ - N. 

D 

[3.02 X 

‘ 

lO”^ + 3.92 X lO"^ + N.(1.72 x lO"^ 


- 6.58 X lO"^ )] , 

(3.19) 





wh«r* T im in °C nnd A « (25 - T). 


The following eapirical relations given by 
Frankenstein and Garner (1967) can be used to coapute the voluae 
fraction of brine in Saline ice - 


10 ^ S ,(- 
si' 


10-^ S,j(. 


10-* S.j(- 


52.56 


45.917 


43.795 


- 2.28) 


+ 0.930) , 


+ 1.189) , 


-0.5 2: T i -2.06 C 

(3.20) 

-2.06 it T 5t -8.2 °C 

(3.21) 

-8.2 it T i - 22.9 °C 

(3.22) 


where is voluae fraction of liquid brine in saline ice and 
is the Salinity of sea-ice aixture. 

The real and iaaginary parts of the coaplex dielectric 
constant of saline ice can be coaputed by the following 
expressions - 



'i 


JB 

, and 

(2.23) 


1 - 3 V. 

D 


SB 

% 'b • 

(3.24) 

In 

the above equations. 

i, si and b refer to ice 


(pure), saline ice and brine. 


3. 3. 3 Dielectric Constant of Snow 

Snow can be divided into two categories : 

(a) dry snow, which is a aixture of ice and air and contains 
no free water, and 

(b) wet-snow, which does contain free water. 
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3. 3. 3. 1 Dry Snow 



3 ^ * * 

where i« the snow density (gn/ca ), and ^ are the real and 
iaaginary parts of the dielectric constant of pure ice and is 
voluae fraction of ice ( « p^ / 0.916 ). The subscripts ds and i 
refer to dry snow and ice respectively. 

The real and iaaginary parts of the dielectric 
constant of dry snow have been calculated using equations (3.25) 
and (3.26) in the frequency range of 1-37 GHz for varying snow 
density and surface teaperature, which is used in the aicrowave 
response study of snow covered terrains. 

The variation of real and iaaginary parts of complex 
dielectric constant with snow density for varying frequencies is 
shown in Figure 3.2. The surface teaperature is taken to be -1 ^C. 

r 

The real part, i.e., dielectric permittivity shows almost 

linear increase with the increase of snow density and is 
independent of frequency. The value of real part of dielectric 

3 

constant increases from 1.2 at snow density of 0.1 gm/cm to 3.1 

3 

when snow density is 0.9 gm/ca . The imaginary part of dielectric 
constant, i.e., dielectric loss factor increases with the 
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Dielectric Loss Factor of Dry Snow (log 



Snow Density (gm/cm^} 

Fi£. 3.2. Variation of real and imaginary parts of complex 

f f I 

dielectric constant of dry snow '^ds with 

with snow density. 



increaa# in frequancy. 


Fifiur* 3.3 ahovs the variation of iaia^inary part of 
dielectric constant with frequency for varying snow density. The 
rate of increase of dielectric loss factor becoaes higher with the 
increase of snow density. 

The variation of loss tangent (tan 6 ) with snow 

s 

density for varying surface temperature is shown in Figure 3.4. 
It is obvious from the figure that the loss tangent of dry snow 
increases with snow density as well as teaperature. The variation 
of loss tangent for different temperatures is smaller at lower snow 
density and becomes significant at higher snow density. 

In summary, the real part of dielectric constant of 
dry snow depends only upon snow density, whereas the imaginary 
part of dielectric constant is highly dependent on snow density, 
frequency, and surface temperature. The snow density 
is directly related to the age of the snowpack, therefore, 
the dielectric constant (both real and imaginary parts) of dry 
snow increases as ageing process prolongs. 


3. 3. 3. 2 Wat, Show 


Modlflmd Dttbye-likm modal has been used to calculate 
the real and imaginary parts of the dielectric constant of wet snow 
and is given by (Hallikainen et al . , 1986)* 


f 

va 




9 9 


A + , and 

1 + ( f / fp )^ 

c ^ / ^0 ) -V* 

1 + ( f / fjj )^ 


(3.27) 


(3.28) 
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wh«r« Im th* relaxation frequency and ia liquid water content 
of wet snow. The constants A, B, C, f^ and x are deternined by 
the following espirical relations - 


A 

m 

1.0 + 1.83 p + 0.02 A. a^°^^ + B. , 


(3.29) 

B 

m 

0.073 A^ , 


(3.30) 

C 

s 

0-073 A 2 , 


(3.31) 

X 

m 

1.31 , 


(3.32) 

^0 

m 

9.07 GHz , 


(3.33) 


ss 

0.78 + 0.03 f - 0.58 x lo"^ f^ , 


(3.34) 

^2 

M 

0.97 - 0.39 f X 10 “^ + 0.39 x lO”^ f^ , 


(3.35) 


s 

0.31 - 0.05 f + 0.87 X 10 “^ f^ , 


(3.36) 

where f 

ia 

the frequency (in GHz) of the signal used. 

^s 

is the 

3 

snow density (ga/ca ). The above expressions are valid 

for 

ranges 

3 :£ f £ 

37 

3 

GHz, 0.09 i 0.38 ga/ca , and 1 S a^ i 

12 % 

. The 


subscript ws refers to wet snow. 

The real and isaginary parts of the dielectric 
constant of wet snow have been conputed and used in all nunerical 
calculations by using above relations for various liquid water 
content, snow density and frequency. 

The behaviour of real and iaaginary parts of 

f ft 

dielectric constant of wet snow and-^ . respectively) with 

frequency is shown in Figure 3.S for various liquid water content 
in the -frequency range of 3*-37 GHz. The density of wet snow is 

3 

taken to be 0.25 ga/ca . The real part of dielectric constant Cor 
dielectric peraittivity) shows a decreasing trend with the increase 
in frequency up to 25 GHz, and afterwards attains a constant value. 
It also increases with the increase in liquid water content of snow 
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•.nd this incr««s« is hifiher at lowsr frequencies. The iaafiinary 
part of dielectric constant (or dielectric loss factor) first 
increases up to about 9 GHz and decreases thereafter continuously 
with the increase in frequency. The rate of increase of ina^lnary 
part of dielectric constant increases with the increase in liquid 
water content of snow. The values of real part of dielectric 
constant of wet snow at 3 GHz for 2% and 10% liquid water content 
are 1.8 and 2.26, and at 37 GHz for the sane liquid water content 
are 1.16 and 1.43 respectively. The values of iaaginary part of 
dielectric constant of wet snow at 3 GHz for for 2% and 10% liquid 
water content are 0.05 and 0.43, and at 37 GHz for the saae liquid 
water content are 0.06 and 0.47 respectively. 

In Fi£ure 3.6, the variation of real part of 
dielectric constant of wet snow with snow density is shown at 9 GHz 
and for different liquid water content. The real part of dielectric 
constant increases linearly with the increase in density and liquid 
water content of snow. The dielectric pernittivity of wet snow 
varies fron 1.18 to 2.7 as the density of snow changes fron 0.10 

3 

to 0.39 gai/coi . 

The real and iaaginary parts of dielectric constants 
of wet snow are highly dependent upon density of snow, liquid water 
content, and frequency. Moreover, liquid water content of snow is 
indirectly related to the surface temperature, therefore, the real 
and imaginary parts of the dielectric constant increase with the 
increase in surface temperature. 

3.3.4. Dl*l«ctrlc Constant of Pure Water 

The real and imaginary parts of dielectric constant of 
pure water has been calculated by using the well known Debye 
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equation (Ulaby «t al . , 1986) and la axprasaed aa - 


at ^ £ 

WO VKO 


V 


1 + (2n£T )‘ 
w 

2n£T ( rt ~ ) 

w '• wO \Kt> ' 


1 + iZnir^y 


where 


^ atatlc dielectric conatant o£ pure water, 
» high-£requency Unit ot a ( ■ 4.9), 

ww V 

s relaxation tine o£ pure water, and 
£ ■< electronagnetic £requency. 


(3.37) 


(3.38) 


The atatic dielectric conatant o£ pure water (-^^q ) 

and the relaxation tine o£ pure water (t ), are £unctiona o£ water 

w 

tenperature and are given by 

2nT (T) » 1.1109 X 10"^° - 3.824 x lO"^^ T + 6.938 

W 

X 10“^^ - 5.096 X 10“^^ , (3.39) 

*wo(^> “ 88.045 - 0.4147 T + 6.295 x lO"* 

+ 1.075 X lO"^ , (3.40) 

where T ia in ®C. The aubacript w re£era to pure water. 

3. 3. S Oinlactrlc Conslamt of Saline Water 

Saline water ia a water containing diaaolved 
aalta. The real and inaginary parta of the dielectric conatant 
of aaline water can be expreaaed aa- 
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mwO awQo 


aw 


^ + 
awca 


, and 


C3.41) 


1 - C2n£T.^)‘ 


2nir ~ « ) 

•w ' «w0 mvco' 


aw 


(3.42) 


1 ♦ 2n«jJ 

whera the aubscript aw refera to aaline water, o*^ ia the ionic 
conductivity of the aqueoua aaline aolution and ® 4.9. 

The dependence of atatic dielectric conatant water 

aalinity S (4») and aurface teaiperature (T^^C) ia given by 

0 w 


* «(T,S ) 

awO'- ’aw-' 


-^awO 


CT,S ) 

V » 0V' 


(3.43) 


where 

*8w0^^’°^ “ 87.134 - 1.949 xio"^ T - 1.276 x 10“^ 

+ 2.491 X lO"* , (3.44) 

and 

a (T,S ) « 1.0 + 1.613 X lO"^ T S_, - 3.656 x lO"^ S„., 

' aw aw sw 

+ 3.210 X lO"® - 4.232 x lO"^ . (3.45) 


The relaxation tine ia given by 

whar* T^.^(T,0) * T CT) («•• aQn. 3.39), and 
aw V 

b (T,S^y) » 1.0 + 2.282 X 10^ T - 7.638 x 10* 

- 7.76 X 10"^ + 1.105 X 10"® S® . (3.47) 

aw aw 

The ionic conductivity for aaline water (<yj^) given by 
Stogryn (1971) takea the for* 
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&AT,S ) 
' av' 


-iC25,S^^) . 


(3.48) 


- 4 , 


f 


wher« *^£(25 , la the ionic conductivity of »ea water at 25 ® C 
and ia given by 

tyiCaS.S^^) « [ 0.18252 - 1.4619 x lo"^ - 2.093 

X 10 "® sl^ - 1.282 X 10 "’ ] . (3.49) 

The function 4> dependa on S and ^ = 25 - T, 

10 

4> * [ 2.033 X 10 “^ + 1.266 x 10“*Zl+ 2.464 

X 10 "^ ,4^ - S^.. ( 1.849 X 10“® - 2.551 

0W 

X lO"’'^ + 2.551 X 10“® )] (3.50) 


The coabination of all the above equationa are valid 
for the aalinity range of 4 to 35 %. . 


3.3.6 Dl«l*ctrlc constant of barer ground 

The dielectric conatant of bare ground been 

computed by using a aeai-eapirical aodel (Oobaon et al . , 1985), 

which ia given by 


a 

& 


1 + 



1) + nl^ (^" - 1) 
' V ' w ' 


(3.51) 


3 . 

where p. ia the bulk soil density (ga/ca ), p ia the density of 

O 00 

3 

solid soil aaterial (» 2.65 ga/ca ) , * ia the coaplex dielectric 

S o 

constant of solid soil aaterial (= 4.7 - j 0), is the liquid 

water content of soil, and « ia the coaplex dielectric constant of 

w 

water. The values of ot and ft are taken to be 0.65 and 1.0 (optiaua 
for all soil types). 
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CHAPTER 4 


RESULTS AND DISCUSSIONS 


4.1 INTRODIICTICM4 

In order to study the potentiality of nicrowave 
radioneters on board and to aid the analysis and interpretation of 
Microwave renote sens ins data, the knowledge of snow and ground 
paraneters affecting the nicrowave renote sensing response is very 
nuch essential. In the present work, the theoretical response of 
nicrowave renote sensing over snow covered terrains is studied. 
The study has been carried out over various two-layered and 
nult i-layered nodels representative of snow covered 

terrains. The results of the detailed nunerical calculations are 
presented to study the potentiality of passive nicrowave radionetry 
to retrieve snowpack paraneters, to nap the snow cover, to nonitor 
snownelt, to estinate the thickness of natural lake ice, and to 
suggest the optinun sensor paraneters to be used in the polar and 
nid- to high-lattitude regions. The variation of nodel paraneters 
of snow as well as the penetration depth o.f microwaves in the snow 
nediun with varying physical paraneters in the frequency range of 
1-37 GHz are also studied. In all the calculations, the real part 
of relative dielectric constant of pure ice is taken 3.15 (as 
mentioned in chapter 3) and imaginary part is taken from Figure 
3.1. The real and imaginary parts of dielectric constants of dry 
and wet snow are taken as discussed in sections 3.3. 3.1 and 


3. 3. 3. 2. 



4. a PARAMErrSH^ GF 

Th« »od«l pi&ra£ft«t <ars of snow, i.«., scattaring (k ), 

£1 

absorption extinction coefficients govern the 

emission behaviour of the snowpack. The snow particles have been 
taken as spherical and uniformly distributed. The oiodel paraneters 
have been computed by utilising Rayleigh expressions for given 
physical properties of the snow sediua like snow density, radius of 
snow particles, liquid water content etc. 

Using equations (2.d) to (2.24), we have studied 

the propagation characteristics of snow. The results are shown in 

Figure® 4. 1-4. 3. Figure 4.1 deasonatrat es the variation of 

scattering CR„)i absorption (k^), and extinction (k *k +k ) 
m e a a 

cctoff icients of dry snow with frequency for snow particles of 1.0 
mm sise. It is evident fro« the Figure 4.1 that the loss of 
eleetrostagnetic energy due to absorption is aore than the loss due 
to scattering up to 19 GHz, and afterwards scattering becomes the 
dominant lose mechanism. The behaviour of model parameters of dry 
snow with the snow particles radius for 10 and 19 GHz is shown in 
Figure 4.2. It is clear from the Figure 4.2 that the scattering, 
absorption and extinction losses are more for 19 GHz than that at 
10 GHz. It can ba seen that the absortion coefficient is 
independent of the snow particles size, whereas scattering 
coefficient is highly dependant on snow particles size. For smaller 
particles, absorption plays an isaportant role, whereas for larger 
pjift i cl iss, scattering becoaer thr primary loss process. 

Thtft variation of absorption, scattering, and 

extinction coef fid eats of wet snow with the radius of snow 
psrtides is shown in Figure 4.3 for 10 and 19 GHz frequency. 
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The liquid v«ter content of enow is taken 1.0% . The behaviour of 
scattering and absorption coefficients is saae as shown in Figure 
4.2. It is clear froa Figure 4.3 that the absorption loss doainates 
over scattering loss in the presence of even a saall aaount of 
liquid water. 

The absorption and scattering losses are highly 
dependent on aicrowave frequencies, and particles size and liquid 
water content of snow. At low frequencies, when the size of snow 
particles is auch saaller than wavelength, absorption plays an 
iaportant role, whereas at higher frequencies when the size of the 
srfow particles is of the order of wavelengths scattering becoaes 
significant. The absorption is the priaary loss process in wet 
snow. Moreover, The total loss of electroaagnetic energy is auch 
higher for wet snow as coapared to that of dry snow. 

4.3 PENETRATION DEPTH OP MICROirAVES IN SNOW MEDIUM 

The penetration depth of aicrowaves in the snow aediua 
is the function of the frequency of aicrowave signal, and particles 
size and liquid water content of snow. In Figure 4.4, the 
variation of penetration depth with density of dry snow (i^=0.0%) 
for various snow particle sizes is shown at 19 GHz. The penetration 
depth of aicrowave signal initially decreases rapidly with the 
increase in snow density and attains a constant depth value, which 
is an optiaua penetration depth. The penetration depths of 

3 

aicrowave signal (F»19 GHz) in dry snow of density 0.1 gn/ca and 
snow particles of radius O.S, 1.0, 1.5 and 2.0 aa are 28.0, 21.0, 
12.3, and 7.0 aeters respectively. 

The penetration depth of aicrowave signal (F»19 GHz) 





with v«t anov danaity ia shown in Figure 4.5. Tha panatration depth 
la atrongly raducad whan tha anew ia vet. The decrease in 
penetration depth la because of tha fact that total loss of 
electronaenetic enarfiy due to scattering and absorption is very 
high in case of wet anov as conpared to that of dry anov. 
Hera, it can be seen that penetration depth renaina sane for 
different particle aizaa up to snow density of 0.19 gn/cn . The 
penetration depth of nicrowave signal at 19 GHz in vet snow is 

3 

approximately 48 ca for anov density of 0.10 gn/cm . 

Figure 4.4 shows the variation of penetration depth in 

3 

vet snow (density » 0.25 gm/cm , and radius of snow particles - 0.5 
ma) with moisture content for varying frequencies. Because of 
the increase in dielectric constant of snow with moisture content 
and consequently total loss of electromagnetic energy, penetration 
depth reduces continuously with the increase in liquid water 
content of snow. The penetration depth also decreases as the 
frequency of microwave signal increases. 

4.4 PASSIVE MICROWAVE REMOTE SENSIM6 RESPf»ISE OF SNCMT 
COVERED TBRRAXMS 

Snow covered terrains can be represented as 
two-layered and multi-layered models.’ The physical and 
electromagnetic properties of each snow layer and underlying ground 
control the microwave remote sensing response. As discussed in 
Chapter 2 , in passive microwave remote sensing, brightness 
temperature is measured from which snow layer characteristics are 
deduced. The detailed study of brightness temperature behaviour has 
been carried out over various two-layered and multi-layered models 
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r«pr*<i«nt«.tiv« of anow cov«r«d terrains in mid- to high-lattitude 
and aountaineous r«||ions. 

■d. 4. 1 Two - Lay«r Models 

A anow covered region can be repreaented as a simple 

‘I 

tvp-layered model. In order to atudy the difference in microwave 
remote aenaing reaponae over bare ground and anow covered ground, 
brightneaa temperature haa been computed with varioua aenaor and 
layered (anow) paraaetera. In Figure 4.7, frequency reaponae of 
brightneaa temperature in vertical and horizontal polarizations for 
microwave aignala incident at 30^ is shown for bare ground and 
ground covered with dry and wet anow. The denaity and average snow 

3 

particles size are 0.25 ga/cm and 1.6 on respectively. The 
brightneaa temperature of ground covered with dry snow shows 
decresing trend with the increase in frequency. For snow-free and 
vet anow covered ground, brightneaa temperature shows increasing 
trend. The brightneaa temperature contrast in vertical and 
horizontal polarizations ia significant at lower frequencies and 
becomes negligible at higher frequencies for wet anow. Horeover, at 
higher frequencies, contrast between the brightness temperature of 
dry anow and wet anow and that of dry anow and bare ground is 
very pronounced. However, it is difficult to separate the areas 
covered with vet snow from snov-free areas in both the 
polarizations. 

The extinction loss of dry anow ia very small 
up to about 19 GHz and microwaves are capable to penetrate through 
the snow layer. Therefore, at low frequencies, brightness 
temperature behaviour ia mainly determined by the underlying ground 
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Fr«qu«ncy (GHz) 


Fig. 4.7. Frequency response of brightness temperature of snow 
covered and bare ground (m^ of dry snow = 0.0%, m^ of wet snow 

= 1.0%, of bare ground = 5%, temperature of dry snowpack = 

272.0 K, temperature of wet snowpack = 275.0 °K, temperature of 

bare ground = 298.0 °K, bulk soil density = 1.7 gm/cm\ 
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and at hifth^r fr«qu«ncl*a by the snow layer. However, in the caae 
of wet enow, brightneee tenperature behaviour ie aainly governed by 
the enow layer. Therefore, a eudden decreaee in brightneee 
tenperature ie encountered when going fron bare ground and ground 
covered with wet enow to the ground covered with dry enow in the 
frequency range of 19-37 GHz. 

The brightneee tenperature contraet in two 

polarizatione (Tgy - Tg^) with frequency ie ehown in Figure 4.8 
for the eane nodele. The brightneee tenperature contraet (Tg^ - 
Tgg) for dry enow ie alnoet frequency independent up to 15 GHz and 
afterwarda it ehowe a decreasing trend. For wet enow, it shows a 
decreasing trend up to 27 GHz before attaining a constant value, 
whereas it decreases slowly in the case of bare ground for the 
entire frequency range. The brightness tenperature contrast (Tgy - 
Tgg) is very snail for wet enow as conpared to that of bare ground. 

Figure 4.9 denonetrates the angular dependence of 
brightneee tenperature at 19 GHz of bare ground and ground covered 
with dry and wet snow for the nodels discussed earlier. The 
brightneee tenperature decreases in horizontal polarization and 
increases in vertical polarization as the incidence angle 
increases. For bare ground, the brightness tenperature in 
horizontal as well as in vertical polarization is highly 
dependent on incidence angle; and for wet snow, the brightness 
tenperature contrast in both the polarizations is insignificant 
even at very high incidence angle. Therefore, to diecrininate the 
snow-free areas fron areas covered with wet snow, incidence angle 
of 50® nay prove to be ideal. 

The negative spectral behaviour of brightness 
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Fig. 4.8. Variation of brightness temperature contrast in vertical 

and horizontal polarizations CT_„ - T„„) with frequency. 

B V JSH 



Fig. 4.9. Angular dependence of brightness temperature for snow 
covered and bare ground. 



t*flip«rature in dry snow im th* unlqua property for diatinguiahing 
it froa wat anew and bar* ground. The polarization infornation way 
iwprove the aeparability between the bare ground and ground covered 
with wet anow. Hence, wapping of anow cover aa well aa 
claaaif ication of anow typea, i.e., the diacriwination of dry and 
wet anow ia poaaible on the baaia of brightneaa temperature 
behaviour ahown in the broad band frequenciea. 

Snowpack uaually undergo ea aeveral welting and 
refreezing cyclea, aa a reault, the aize of anow particlea geta 
changed throughout ita life hiatory. The aize of anow particlea 
increaaea aignif icantly with alight wetting. In Figure 4.10, the 
frequency reaponae of brightneaa tewperature for varying anow 
particle aizea and anow layer thickneaaea ia ahown at 5(f look 
angle and horizontal polarization. The brightneaa tewperature 
rewaina conatant up to about 9 GHz for anow particlea of 0.8 and 
1.2 aw radii, and decreaaea aharply thereafter before attaining 
a conatant value. Aa ahown in Figure 4.10, the brightneaa 
tewperature decreaaea aa the particle aize and depth of anow layer 
increaae. For very awall particlea, brightneaa temperature ia leaa 
senaitive to anow thickneaa. 

The brightneaa temperature behaviour with anow 
particlea aize ia ahown in Figure 4.11 for 15, 19, 25 and 37 GHz 
and at look angle of 50^. Aa diacuaed in the aection 4.2, 
acattering loaaea dominate a dry anow cover aa the frequency and 
anow particlea aize increaae, which cauaea the reduction in 
brightneaa tewperature. The decreaaing trend of brightneaa 

temperature increaea with the increaae of frequency. For the anow 
layera containing particlea of 0.2 and 2.0 wa radii, the brightneaa 
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tdaperaturea are approxlnately 254.0 and 180. O^K reapectively at 15 
GHz and are 254.0 and 5.0 reapectively at 37 GHz. Fron Figure 
4.11, it ia concluded that for very aaall anov particlea, 
the brightneaa teaiperature ia independent of frequency, whereaa for 
larger particlea it ahowa very atrong frequency dependence. 

To eatiaate the daily anov cover extent and accurate 
runoff froa a anov covered area, knowledge of anov thickneaa ia 
very iaportant (Schultz, 1988). Figure 4.12 ahowa the behaviour of 
brightneaa teaperature with anow thickneaa for 19, 25, and 37 GHz 

and look anglea of 30^ and 50^. The average anow particlea aize ia 
taken 1.0 aa. At 19 GHz, brightneaa teaperature ia alaoat 
independent of the thickneaa of the anow layer and aa the frequency 
increaaea, the negative alope of the curve increaaea. The 
brightneaa teaperature decreaaea aa the look angle increaaea froa 
30^ to 50^. However, the behaviour of brightneaa teaeprature ia 
aiailar for 30^ and 50*^ look anglea. It ia evident froa the Figure 
that for anov layer containing particlea of 1 aa aize, aore 
accurate eatiaation of anov thickneaa can be aade at 37 GHz. 

The effect of anov particlea aize and thickneaa on 
brightneaa teaperature at 37 GHz and 3(f* look angle ia ahown in 
Figure 4.13. The brightneaa teaperature firat decreaaea with the 
increaae in anov thickneaa and becoaea conatant after a aaturation 
depth for a particular aize of anov particlea. For larger anow 
particlea, acattering loaa of electroaagnetic energy becoaea ao 
aignif leant that even a aaall increaae in the thickneaa of anov 
layer reaulta in a draatic decreaae in brightneaa teaperature. 
Therefore, thickneaa of anow containing larger particlea can not be 
Inferred at 37 GHz. For anow particlea of 1.0 aa aize, the 
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snow layer thickness upto about 2.5 m can be estimated accurately. 

Therefore, the thickness of snow layer and particle 

size strongly influence the emission behaviour of the snowpack. The 

determination of snow water equivalent, i.e., how much water can be 

d 

derived from a snow covered area I- X p dz, where d is the 

0 ® 

thickness of the snow layer) is very important for the prediction 
of runoff from a river basin. The penetration depth of 
microwaves in wet snow is very small (as discussed in section 
4.3), therefore, the thickness of wet snow can not be estimated. 
The size of snow particles generally varies between 1 to 2 mm. 
Therefore, for average size of snow particles, frequency range of 
25—37 GHz will be suitable for the extraction of snow thickness 
information. 

The snow wetness parameter is very important in 
estimating water reserve for irrigation and hydropower (Schultz, 
1988). The variation of brightness temperature as a function of 
frequency for various snow water content is shown in Figure 4.14 in 
horizontal and vertical polarizations. For vertical polarization, 
brightness temperature remains nearly constant up to approximately 
30 GHz and is followed by a slight decrease. In case of horizontal 
polarization, the brightness temperature increases sharply up to 27 
GHz and decreases thereafter. The decrease in brightness 
temperature may be attributed to the dominating behaviour of 
scattering as loss mechanism at higher frequencies. The variation 
in brightness temperature due to the change in liquid water content 
of snow is hardly discernable in vertical polarization. Therefore, 
vertical polarization is not at all suited for the estimation of 
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Fifi. 4.14. Frequency response of brightness temperature of wet 
snow with liquid water content of snow as a parameter. 



Fig. 4.15. Brightness temperature versus liquid water content of 
snow for different frequencies. 
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water content of snow. However, the brightness temperature 
variation for different water content is significant in horizontal 
polarization. 

In order to estimate the free liquid water content in 
snow, the effect of snow water content on brightness temperature is 
shown in Figure 4.15 for various frequencies. The average size of 
the snow particles is taken to be 1.0 mm. It is obvious from the 
Figure that the brightness temperature of dry snow remains lower 
than that of wet snow up to about SX liquid water content (except 3 
(3Hz> and is higher for further increase in snow wetness. For 3 GHz, 
the brightness temperature initially increases upto 2.5% snow water 
content and then decreases continuously, whereas for higher 
frequencies brightness temperature increases only upto 1% snow 
water content and decreases afterwards. It is evident from the 
Figure 4.15 that a particular brightness temperature represents two 
values of snow water content. Koreover, at lower frequencies the 
brightness temperature of dry snow becomes higher than that of wet 
snow with higher liquid water content. To resolve this ambiguity, 
multi-frequency approach should be preferred. At higher frequency, 
the decreasing trend of brightness temperature with the increase in 
liquid water content (after 1%) decreases. Therefore, in order to 
estimate the liquid water content of snow precisely, lower 
frequencies must be used. 

The information of onset of snow melt can also be 
made available, which is a very important parameter for the 
prediction of timing of runoff from a river basin. Snow becomes wet 
during the day time and becomes dry in the night time. This 
phenomenon can be used to detect the onset of snow melt because a 
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large diurnal fluctuation in brightness tenperature will be 
observed at higher frequencies, especially at 37 GHz (froa 
Figures 4.7 and 4.14). 

Generally, the interfaces between air and snow aediun, 
and snow and ground aediua are seldom saooth due to the presence of 
snow melting and local winds. In order to study the effect of 
surface roughness on microwave remote sensing response, brightness 
temperature as a function of frequency is shown in Figure 4.16 for 
smooth and rough (<? s i cm) surfaces of dry and wet snowpacks. The 
average temperatures of dry and wet snowpacks are taken to be 372 
and 275 respectively. The brightness tenperature is higher 
when the upper surface of dry snow is rough in 3-18 GHz frequency 


range and 

is 

lower after 

18 GHz. In 

the 

case of wet 

snow, the 

variation 

in 

brightness 

temperature 

is 

signi f icant 

at lower 


frequencies and is hardly discernable at higher frequencies. As 
shown in Figure 4.16, the surface roughness effect in case of dry 
snow is minimised at 18 GHz. 

The angular dependence of brightness temperature at 10 
GHz is shown in Figure 4.17 for the same model parameters as 
discussed in Figure 4.16. The variation in brightness tempearature 
in horizontal polarization for plane and rough surfaces increases 
as the incidence angle increases, whereas in the case of vertical 
polarization it shows just opposite behaviour to that of horizontal 
polarization. It is clear from the Figure 4.17 that the effect of 
surface roughness is almost eliminated in vertical polarization at 
about 50^ look angle. Koreover, brightness temperature contrast 
(Tgy - Tgg) decreases as roughness is introduced. 

To study the roughness effect of top and bottom 
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Frequency (GHz) ^ 

Fig. 4.16. Spectral variation of brightness temperature for smooth 
and rough top snow surfaces. 
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Incidence angle (degrees) 


Fig. 4.17. Angular dependence of brightness temperature for smooth 
and rough top snow surfaces. 







surfaces of dry snowpack, angular dependence of emissivity is shown 
in Figure 4.18 at 5 Qlz frequency. It shows that at lower frequency 
<i.e., when penetration depth is large) bottom surface topography 
plays an important role in governing the emission behaviour of the 
snow packs. 

4.4.8 MULTI-LAYER MOOBLS 

Variable exchange of heat and mass on the snow surface 
leads to accumulation of snow in layers. The physical 
characteristi cs of these layers differ from one to another. As the 
thickness of snow increases, the bottom layers become more compact 
and dense. Due to equilibrium thermodynamics, crystals of small 
sizes and of spherical shape are developed and when there is a 
steep temperature gradient with depth, large crystals are 
developed. Here, the passive microwave remote sensing response of 
multi-layered models representative of mountaineous and polar 
regions is discussed for various physical conditions. In all the 
models, the temperatures of dry and wet snow layers are assumed to 
be S72.0 and S75.0**K respectively. 

Figure 4.19 shows the frequency response of brightness 
temperature in horizontal and vertical polarizations for the models 
given in the legend. The moisture content and thickness of each 
snow layer are taken to be 0 . 07 , and 10 cm respectively. The 
brightness temperature in both the polarizations remains same for 
the two models up to 11 Qiz and is higher afterwards for the model 
(1) in which the snow particles size decreases with depth as 
compared to the model (2) in which the snow particles size 
increases. The brightness temperature contrast between the two 
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Fig. 4.1S Effect of top and bottom surface roughness of snow on 
emisalvity . 
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Fig. 4.19 Spectral variation of brightness temperature for multi- 
layered models with decreasing and increasing snow particles 
size with depth. 
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models is higher in vertical polarization than horizontal 

polarization. 

The frequency response of brightness temperature in 
vertical and horizontal polarizations over models representative of 
yearly and multiyear snowpacks is shown in Figure 4.20. The 

yearly and multiyear snow layers are charactar ized with varying 
density and particle size of snow. The moisture content and 

thickness of each snow layer are taken to be the 0.0% and 20 cm 
respectively. The behaviour of brightness temperature is similar as 
earlier. In horizontal polarization, as the number of layers of the 
snowpacks increases, the brightness temperature decreases 
significantly for the entire frequency range. However, in vertical 

polarization, the brightness temperature remains same for the 

models (1) and (2) up to about 13 GHz frequency and as the 
frequency increases, snowpacks containing lesser number of layers 
show higher brightness temperature. The brightness temperature 
contrast for the three models increases with the increase of 

frequency in vertical polarization. 

The angular dependence of brightness temperature for 
the same models (as in Figure 4.20> is shown in Figure 4.21. The 
brightness temperature decreases with the increase of incidence 
angle in horizontal polarization, whereas it increases in vertical 
polarization. The brightness temperature contrast in both the 

polarizations (Tgy - Tgj^J decreases with the decrease in the number 
of snow layers, i.e., multiyear snowpacks show higher brightness 
temperature contrast. Therefore, it can be concluded that the use 
of proper frequency and polarization measurements can reveal 
information about the age of accumulated snow. 
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Snow undergoes melting and refreezing cycles with time 
in response to change in weather conditions. There are two types of 
melting and refreezing cycles — (i) when the upper layers of the 
snowpack are wet and as the temperature falls, the whole snowpack 
refreezes, and (ii) when the whole snowpack is wet and the upper 
layer refreezes during cold nights. In order to monitor the melting 
and refreezing cycles, frequency response of brightness temperature 
in horizontal polarization is shown in Figure 4.22. The physical 
parameters of each snow layer are given in the legend of the 
Figure. In the melting and refreezing cycle of type 1, i.e., when 

the upper layers of the snowpack are wet, brightness temperature 
increases sharply up to 9 GHz and this rate decreases with further 
increase of frequency in horizontal polarization. But as the 
temperature falls and snowpack refreezes, brightness temperature 
first remains constant up to 15 GHz and decreases rapidly 
afterwards. In the case of refreezing and melting cycle of type 2, 
i.e., when the whole snowpack is wet, brightness temperature 
follows the same behaviour as the snowpack of early stage of type 
1. In the later stage of the cycle of type 2, as the upper layer 
refreezes, brightness temperature first shows increasing trend up 
to 15 GHz and decreases rapidly thereafter. The brightness 
temperature of snowpack containing only one layer of snow with 1 % 
liquid water content is higher (however, variation in brightness 
temperature is small) as compared to that of wet multi-layered 
snowpacks (2 and 3). In vertical polarization, almost same 
behaviour is observed (shown in Figure 4.23). However, in vertical 
polarization, brightness temperature is less sensitive to the 
changes that occur during the two types of melting and refreezing 
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cycles. 


The angular dependence of brightness temperature is 
shown in Figure 4-24 at 19 GHz frequency for the models given in 
the legend. The multi-layered models (1) and (2) are the 
representative of the early stage of melting and refreezing cycles 
of type 1 and S respectively, whereas model (3) is a representative 
of yearly snowpack in which the snow layer is wet. The brightness 
temperature contrast in horizontal and vertical polarizations 
(Tgv “ very small for partially or completely wet snowpacks 

as compared to that of dry snowpacks (from Figures 4.21 and 4.24). 
This difference is almost same for the models 1, 2 and 3 (Figure 
4.24). It can be inferred from the Figure 4.24 that the brewster 
angle (the angle at which the brightness temperature is maximum) is 
lower for multiyear snowpacks as compared to that of yearly 
snowpacks. 

In order to study the effect of roughness of 

multi-layered snow surfaces, the brightness temperature behaviour 
with frequency is shown in Figure 4.25 for the models given in the 
legend at 50° look angle. The surface height variation of rough 
surface is taken to be 1 cm. For the snowpack, in which the top 

surface is rough, brightness temperature is higher as compared to 
that of smooth surface in 1—37 GHz frequency range. The variation 

in brightness temperature due to rough surface is pronounced in 

horizontal polarizat ion, whereas it remains same in vertical 
polarization at 50° look angle. 

The angular dependence of brightness temperature at 25 
GHz frequency for dry and wet snowpacks in both the conditions of 
plane and rough surfaces is shown in Figure 4.26. The variation in 
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Fie- 4.24 Angular variation of brightness temperature of single 
and multi-layered wet snowpacks with varying wetness conditions 
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Fig. 4.26 Angular variation of brightness temperature for plane 
and rough top surfaces of dry and wet multi-layered anowpacks . 
















brightness temperature for smooth and rough surfaces increases with 
the increase of incidence angle for dry snowpacks in horizontal 
polarization and becomes negligible after 40° look angle in 
vertical polarization. In case of wet snowpack, the variation in 
brightness temperature for smooth and rough surfaces increases with 
the increase of incidence angle in both the polarizations. However, 
the polarization difference in brightness temperature ) 

decreases as the surface of the snowpack becomes rough. 

Therefore, the effect of surface roughness on 

microwave remote sensing response can be suppressed by the proper 
selection of sensor parameters such as frequency, incidence angle 
and polarization. 

4,. 5 PASSIVE MICROWAVE REMOTE SEN5IMO RESPCMSE OP LAKE ICE 

In polar and high— latti tude regions, the surface of 
the lake is often covered by an ice layer. The lake water is seldom 
pure as it contains dissolved salts and minerals as well as various 
chemicals, which increases its salinity. In order to estimate the 
thickness of ice layer overlying lake water, variation of 

brightness temperature as the function of ice thickness in 

horizontal polarization is shown in Figure 4.27 at the look angle 

of 50° for varying frequencies in both cases of pure and saline 
water. The average physical temperatures of ice layer and that 
of lake water are taken to be 272°K. The brightness temperature of 
lake water Ci.e. when d=0, where d is the thickness of ice layer) 
is lower than that of lake ice. The brightness temperature 
increases with the increase of frequency and ice thickness. The 
brightness temperature shows increasing trend with increase of 
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Fig. 4.28. Angular response of brightness temperature of lake ice 
at 5 and 10 GHz frequency. 
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Bright. Temperature of Saline Ice (Kelvin) 




frequency, the trend is more pronounced for saline ice than for 
pure ice. The brightness temperature of saline ice is higher as 
compared to that of pure ice for a particular frequency and ice 
thickness. For saline ice, the brightness temperature first 
increases rapidly with the increase of ice thickness and becomes 
constant after an optimum thickness for a particular frquency. 
Usually, lake water is never pure, therfore, it is clear from the 
Figure 4.27 that the frequency of 5 GHz will be most suitable for 
the estimation of lake ice thickness. 

The angular dependence of emissivity of lake ice at 5 
and 10 GHz is shown in Figure 4.28. The emissivity of pure and 
saline ice decreases with the increase of incidence angle in 
horizontal polarization, and increases in vertical polarization. 

The variation in brightness temperature for pure and saline ice 
decreases with the increasing look angle in horizntal polarization. 
Therefore, incidence angle of 0° will prove to be more 

advantageous . 

Generally, the surface of lake ice is overlain by the 
snow layer, which modifies the emission behaviour of lake ice. If 
the snow layer is dry, of low density and has small ice particles, 
it will appear transparent electromagnetically for microwaves. 
However, the presence of wet snow will mask the emission from the 
underlying ice layer. Therefore in sensing the lake ice thickness, 
lower frequencies must be used to suppress the effect of snow layer 
on top. 
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CaUPTER s 


CONCLUSION 


The snowfall is a frequent phenomenon for most of the 
year in the mid— and high-lattitude, and mountaineous regions. Due 
to continuous snowfall, snow accumulates on the surface and- forms a 
thick layer. In the polar regions, the yearly accumulated snow 
piles up and forms a large mass of snow and ice. The physical and 
electromagnetic characteristics of the snow or ice varies with 
depth. Timely information about snow and ice cover is needed to 
predict the water supply, to anticipate flooding, to forecast crop 
yields and estimate freeze damage, to monitor climate, to obtain 
insights advantageous to military and foreign policy, and to guide 
winter sport activities. In this study, theoretical microwave 
remote sensing response of snow covered terrains has been computed. 
The numerical results show that the passive microwave remote 
sensing response of snow covered areas is highly dependent on the 
frequency, incidence angle and polarization of the microwave 

signal, and on the snowpack parameters. 

The knowledge of dielectric constant is very important 
for the proper interpretation of microwave remote sensing data. The 
real and imaginary parts of complex dielectric constant of dry snow 
have been computed by using Polder-Van Santen two— phase mixing 
formula and that of wet snow by modified Debye-like model. The real 
part of dielectric constant of dry snow is the function of snow 
density and increases with the increase in snow density. The 
imaginary part increases with the increase of frequency. 



temperature and snow density. The real and imaginary parts of 
dielectric constant of wet snow are the function of frequency, snow 
density and moisture content of snow. The imaginary part of 
dielectric constant of wet snow is considerably higher as compared 
to that of dry snow. 

Scattering is the primary loss process in dry snow, 
whereas absorption becomes the dominant loss mechanism in wet snow 
due to significant increase in the imaginary part of dielectric 
constant when liquid water is present. The penetration depth of 
microwave signal depends upon frequency, snow density, snow 
particle size and free liquid water content of snow. The 
penetration depth is drastically reduced in the presence of even a 
small amount of liquid water and decreases with the increase in 
frequency, particle size and density of snow. 

The brightness temperature of dry snow shows negative 

spectral slope, whereas wet snow and bare ground show positive 

slope. The wet snow can be discriminated from bare ground by 
measuring the brightness temperature contrast in vertical and 
horizontal polarizations at 50° incidence angle. It is concluded 
that the multifrequency and multiplarization approaches can be used 
for the snow characterization. 

The snow particle size increases significantly as the 
snow layer becomes wet, which creates ambiguity in estimating the 
snow thickness. For the average snow particles size, frequency 
range of 25—37 GHz is most suited for obtaining the dry snow 

thickness information. The thickness of wet snow can not be 
determined due to the limited penetration depth of microwave 

signals in the presence of water. The diurnal fluctuation in 
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brightness temperature can be used to detect the onset of snowmelt, 
which is a very useful parameter for predicting the timing of 
runoff from a river basin. The free liquid water content of snow 
can be estimated accurately at lower frequencies. The effect of 
roughness is almost eliminated in vertical polarization at 50° look 
angle . 

Snow generally accumulates in a layered structure and 
physical parameters of each layer change with depth and time. 
Therefore, the dielectric constant is also changed accordingly, 
which causes the brightness temperature of multi-layered models to 
change when compared to that of two-layered models. By the use of 
proper frequency and incidence angle, the number of layers present 
as well as the characteristics of each layer in a snowpack can be 
determined which can be used in monitoring melting and refreezing 
cycles of snow covered terrains. Our results also show that the 
lake ice thickness can be estimated by using microwave signal of 
about 5 GHz frequency. 

The present study of microwave remote sensing of snow 
covered terrains suggests that the multi-channel microwave 
radiometers operating near 37, S5 and 5 GHz may prove to be ideal. 
The look angle of 50° seems to be reasonable for determining the 
the maximum information about the snowpacks. 

The Indian Space Research Organization, Bangalore is 
planning to launch remote sensing satellites with microwave sensors 
on board in early nineties. The results discussed here will be 
useful in the development of appropriate microwave sensors and in 
the accurate interpretation of microwave remote sensing data of 
snow covered regions. 
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The data on dielectric constant of icer and dry 
and wet snow are almost non-existant in microwave frequency 
range. Therefore, detailed experimental studies must be carried out 
to study the behaviour of dielectric properties of pure ice, and 
dry and wet snow in the VHF and microwave ranges. 
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APPENDIX - A1 


Q* * * ic t * * t in it * t * fc it 1 1 1 * ic t t It t * * t* it it -k * 1 1 fcit * * He t * ie t * -k t * It * "k * * t ie * * ic t * 1 1 

c****** this program computes the real and imaginary parts 

C****** OF DIELECTRIC CONSTANT OF DRY SNOW / PROGRAMMED BY 
C****** SUSHIL K.SRIVASTAV / 

O********** ***************************************************** 

c****** VARIABLES INVOLVED ARE FREQUENCY, REAL AND IMAGINARY 
C****** PARTS OF DIELECTRIC CONSTANT OF PURE ICE, AND DENSITY 
C****** OF DRY SNOW AND PURE ICE. 

Q** ************************** t********* ****************** ******* 

C****** PI » FREQUENCY IN GHz 

C****** EPICPX « DIELECTRIC LOSS FACTOR OF PURE ICE 
C****** RHOS = DENSITY OF DRY SNOW 
C****** RHOI = DENSITY OF PURE ICE 
C****** VI a VOLUME FRACTION OF ICE 

C****** EPIRL = DIELECTRIC PERMITTIVITY OF PURE ICE 
C****** EDSRL a DIELECTRIC PERMITTIVITY OF DRY SNOW 
C****** EDSeX a DIELECTRIC LOSS FACTOR OF DRY SNOW 

Q*************************************************************** 

DIMENSION F1(100),EPICPXC100) 

REAL FI , F , RHOS , RHO I , VI , EDSRL , BC , CC , DC , EC , EDSCX 
OPEN(UNIT=26,FILE=*M5.DAT’ ) 

OPEN(UNITa27,FILE='F4.RES' ) 

PIa3. 142857 

RHOIaO.916 

EPIRL=3.15 

DO 200 IRSal0,50,10 

DO 100 Ial,37 

READC26,*)F1CI),EPICPX(I) 

FaFl(I)*1.0E09 
RHOSa IRS/ 100.0 
VlaRHOS/RHOI 

EDSRLa(1.0+0.508*RH0S)**3 

BC=EDSRL**2 

CCa 2.0* EDSRL 

DC=3.15+CC 

ECa3.15+(2.0*BC) 

EDSCPXa3.0*VI*EPICPX(I)*BC*(CC+1.0)/(DC*EC) 

WR I TE C 2 7 , * ) EDSRL , EDSCPX 
100 CONTINUE 

200 CONTINUE 

STOP 
END 


Q****** 
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APPENDIX > A2 


Q*1e******l!*1c1c*1t*****1t1i*******1c**1t1t**1i***1c****1it*1c***1c******1i1t**** 

C****** the PROGRAM COMPUTES THE REAL AND IMAGINARY PARTS 
C****** OF THE COMPLEX DIELECTRIC CONSTANT OP WET SNOW 
C****** by using DEBYB-LIKE model / PROGRAMED BY SUSHI L K. 

C****** SRIVASTAV / 

C****** VARIABLES INVOLVED ARE FREQUENCY, DIELECTRIC LOSS FACTOR 
C****** OF PURE ICE, SNOW WETNESS, AND DENSITY OF WET SNOW. 

Q*t ************************* **************** t* ************ ******* 

C****** FI =» FREQUENCY IN Hz 

C****** EPICPX = DIELECTRIC LOSS FACTOR OF PURE ICE 
C****** MV = VOLUMETRIC WATER CONTENT IN PERCENTAGE 
C****** RHOS » DENSITY OF WET SNOW IN C.G.S. UNITS 
C****** RHOI = DENSITY OF ICE IN C.G.S. UNITS 
C****** EWSRL » DIELECTRIC PERMITTIVITY OF WET SNOW 
C****** EWSCX = DIELECTRIC LOSS FACTOR OF WET SNOW 

Q************************************* *************************** 

DIMENSION F1(200),EPICPX(200),EWSRL(200),EUSCX(200) 

REAL FI, EPICPX, FREQ, RHOS, RHOI, MV, A1,A2,B1, A, B,C,X,F0,Z, EWSRL, El 
OPEN(UNIT=21,FILE='Al.DAT’) 

OPEN(UNIT=25,FILE='A6.RES’ ) 

PI=3. 142857 

RHOS=0.25 

RH0I=0.916 

DO 400 IM=1,9,2 

DO 300 I»l,37 

READ(21,*)F1CI),EPICPXCI) 

MV=IM/1.0 

FREQ=37.0 

A1=0. 78+0. 03*FREQ-0- 58*1, 0E-03*(FREQ**2) 

A2=0. 97-0. 39*FREQ*1.0E-02+0. 39*1 .OE-03* (FREQ* *2) 

B1=0. 31-0. 05*FREQ+0. 87*1. 0E-03*(FREQ**2) 

A=1.0+(1.83*RHOS)+C0.02*A1*(MV**1 .015))+B1 

B=0.073*A1 

C=0.073*A2 

X=1.31 

F0*9.07 

Z=FREQ/F0 

EWSRL(I)=A+C(B*(MV**X))/(1.0+Z*Z)) 

EWSCX(I)=(C*Z*MV**X)/(1.0+Z*Z) 

WRITE(25,*)IM,EWSRL(I),EWSCX(I) 

300 CONTINUE 

400 CONTINUE 

STOP 
END 
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APPENDIX - A3 


Q ******* ****11 »**** ****************** ************** *************** 

c****** ihE program computes the real and imaginary part 

c****** OF THE COMPLEX DIELECTRIC CONSTANT OF PURE UATBR 

C****** by using DEBYE EQUATIONS / PROGRAMMED BY SUSHI L K. 

C****** SRIVASTAV / 

***************************************************************** 
C****** VARIABLES INVOLVED ARE FREQUENCY AND TEMPERATURE. 

Q** ******************************************************** ****** 

C****** Fl= FREQUENCY IN GHz 

C****** T= PHYSICAL TEMPERATURE OF WATER 

C****** EPWRL = DIELECTRIC PERMITTIVITY OF PURE WATER 

C****** EPWCX = DIELECTRIC LOSS FACTOR OF PURE WATER 

Q ************************* *************************************** ******1 

DIMENSION F1C100),EPWRL(100),EPWCX(100) 

REAL F1,F,T,EWINF,EWO,X,A,B,TOW,C,D,EPWRL,EPWCX 
OPEN(UNIT=22, FILE=’ FREQ. DAT’ ) 

OPENCUNIT=23,FILE=’A2.dat' ) 

PI=3. 142857 
DO 100 IT»1,26 
DO 100 1=1,28 
T=IT-1.0 
READC22,*)F1(I) 

F=F1(I)*1.0E09 

EWINF=4.9 

EW0=88.045-0.4147*T+6.295E-04*T*T+1.075E-5*T**3 

X=2.0*PI 

A=1 .1109E-10~3.824E-12*T 

B=6.938E-14*T*T-5.096E-16*T**3 

TOW=(A+B)/X 

C=(X»F*TOW)**2 

D=1.0+C 

EPWRLC I )=EWINF+ ( (EWO-EWINF) /D) 

EPWCX(I)=(X*F*TOW*(EWO-EWINF))/D 

WRITE(23,*)F1CI),EPWRL(I),EPWCX(I) 

100 CONTINUE 

200 CONTINUE 

STOP 
END 
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£**** ******* It******** ******************************************** 

c****** this program computes the real and imaginary part 

c****** OF THE COMPLEX DIELECTRIC CONSTANT OF SALINE 
c****** water by using STOGRYN'S FORMULATION / PROGRAMMED 
C****** BY SUSHI L K. SRIVASTAV / 

Q**************************************************************** 

c****** VARIABLES INVOLVED ARE FREQUENCY, TEMP. AND SALINITY. 

Q******************************* ********************************* 

c****** Fl» FREQUENCY IN GHz 

C****** BPO a PERMITTIVITY OF FREE SPACE 

C****** T a PHYSICAL TEMPERATURE OF SALINE WATER 

C****** s= SALINITY OF WATER 

C****** EPWRL a PERMITTIVITY OF SALINE WATER 

C****** EPWCX a DIELCTRIC LOSS FACTOR OF SALINE WATER 

Q*** ****************************************************** ******* 


DIMENSION F1(100),EPWRLC100),EPWCXC100) 

REAL T , S , FI , F , EPO , EPTO , A , B , ATS , EPTS , TM , TO , C , D , BTS 
REAL TOTS , DEL , E , G , PHI , SIG , SIGMI . XX , EPWRL , EPWCX 
OPEN(UNIT=24,FILE=' FREQ. DAT’ ) 
OPEN(UNIT=25,FILEa’Bl.DAT’ ) 

PI=3. 142857 
EP0=8.854E“12 
DO 200 1=1,28 
T=25.0 
S=35.0 

READ(24,*)F1(I) 

F=F1(I)*1.0E09 

EPT0=87. 134-1. 949E-1*T-1.276E-2*T*T+2.491E-4*T**3 

Aal.0+1.613E~5*T*S-3.656E-3*S 

B=3.210E-5*S*S-4.232E-7*S**3 

ATSaA+B 

EPTSaEPTO*ATS 

TMal.ll09E-10-3.824E-12*T+6.938E-14*T*T-5.096E-16*T**3 

T0=TM/C2.0*PI) 

Cal.0+2.282E-5*T*S-7.638E-4*S 

Da? . 760E-6*S*S-1 . 105E-8*S**3 

BTSaC-D 

TOTS=TO*BTS 

DELa25.0-T 

Ea2 . 033E~2+1 . 266E-4*DEL+2 . 464E-6*DEL*DEL 

GaS*(1.849E~5-2.551E-7*DEL+2.551E-8*DEL*DEL) 

PHIaDEL*(E-G) 

SIG=S* (0.18252-1. 4619E-3*S+2.093E-5*S*S-1.282E“7*S**3) 
SIGMIaSIG*EXP(-PHI) 

XX=2.0*PI*F*TOTS 

AAl a ( XX* ( EPTS- 4 . 9 ) / ( 1 . 0+XX*XX ) ) 
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200 


BB1=(SIGMI/(2.0*PI*F*EP0)) 

EPURL(I)=4.9+((EPTS-4.9)/(1.0+XX*XX)) 

EPUCXCI)=AA1+BB1 

WRITEC25 , * )F1 ( I) , EPWRL( I ) , EPWCX( I ) 

CONTINUE 

STOP 

END 
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C********* ******************************************************* 

c****** the program computes the absorption, scattering and 
c****** extinction coefficient, and depth of penetration of 
c****** microwaves in dry and wet snow by using RAYLEIGH 

C****** EXPRESSIONS / PROGRAMMED BY SUSHIL K. SRIVASTAV / 

( 2 ** ************************************************** ************ 
c****** VARIABLES INVOLVED ARE FREQUENCY, TEMPERATURE, REAL AND 
C****** IMAGINARY PART OF COMPLEX DIELECTRIC CONSTANT OF PURE ICE, 
C****** DENSITY OF SNOW AND PURE ICE, SNOW PARTICLE SIZE, AND 
C****** SNOW WETNESS. 

Q* t * * 1 1 * ic * 1 1 it it t it it 1 1 1 * 1 1 ic 1 1 ie * * ic ie ic 1 1 * * it 1 1 * * t * * 1 1 ic * ic t * * ic * t t * t * * -k 

C****** FI » FREQUENCY IN Hz 

C****** EPICPX = DIELECTRIC LOSS FACTOR OF PURE ICE 

C****** Rss RADIUS OF SNOW PARTICLES IN METRES 

C****** MV « SNOW WETNESS EXPRESSED IN PERCENTAGE 

C****** RHO:i = DENSITY OF PURE ICE IN C.G.S. UNIT 

C****** RHOS « DENSITY OF SNOW IN C.G.S. UNIT 

C****** VI - VOLUME FRACTION OF ICE 

C****** EDSRL = DIELECTRIC PERMITTIVITY OF DRY SNOW 

C****** EDSCPX = DIELECTRIC LOSS FACTOR OF DRY SNOW 

C****** EWSRL = DIELECTRIC PERMITTIVITY OF WET SNOW 

C****** EWSCPX * DIELECTRIC LOSS FACTOR OF WET SNOW 

C****** eFB = COMPLEX DIELECTRIC CONSTANT OF BACKGROUND MEDIUM 

C****** EPWRL = DIELECTRIC PERMITTIVITY OF PURE WATER 

C****** EPWCX = DIELECTRIC LOSS FACTOR OF PURE WATER 

C****** ePSR = DIELECTRIC CONSTANT OF SNOW RELATIVE TO THAT 

C****** OF THE BACKGROUND MEDIUM 

C****** KAI = ABSORTION COEFFICIENT OF ICE PARTICLES 

C****** KAB = ABSORPTION COEFFICIENT OF BACKGROUND MEDIUM 

C****** KA = TOTAL ABSORPTION COEFFICIENT (KAI+KAB)OF SNOW MEDIUM 

C****** KS» SCATTERING COEFFICIENT OF SNOW MEDIUM 

C****** KE = EXTINCTION COEFFICIENT OF SNOW MEDIUM 

C****** DELTP = PENETRATION DEPTH 

C**********************************************************************i 

DIMENSION F1C500) ,EPICPX(500) 

REAL FI, FREQ,R,Rl,RHOI, RHOS, VI, EDSRL, BC.CC, DC, EC, EPICPX 
REAL EDSCPX,EWSRL,EWSCPX,A2,B1,A,B,C,X,X1,F0,Z,CA,DA 
REAL 2, EWINF.EWO, TOW, EPWRL, EPWCX, ALB, DELTP 
REAL K1,K2,KA,KE,KS,L,NV,LMDA0,KAI,KAB,KA1,KS1,KE1 
COMPLEX EP2,AK,EPSR,EPB,EPW,A1 
OPEN(UNIT=21.FILE= ’SKS1.DAT’) 

0PEN(UNIT=25,FILE=’SKS1 .RES’ ) 

PI=3. 142857 

RHOI=0.916 [ 

DO 500 IM=1,13,1 
DO 400 1=1,37 



PAGE : 7 


READ(21,*)F1(I),EPICPX(I) 

FREQ=F1(I)/1.0E09 

MV=(IM-1.0) 

DO 300 IR=20,500,10 

R=(IR/100.0)*1.0E~03 

R1=R*1.0E03 

DO 200 IRS«10,39,1 

RHOS=(IRS/100.0) 

K1=(2.0*PI*F1(I))/3.0E08 

VI=RHOS/RHOI 

IF(MV.GT.O.O)GO TO 15 

EDSRL=(1.0+0.508*RHOS)**3 

BC=EDSRL**2 

CC=2«0*EDSRL 

DC=3.15+CC 

EC»3.15+(2.0*BC) 

EDSCPX=3.0*VI*EPICPX(I)*BC*(CC+1.0)/(DC*EC) 

EP2=CMPLX(BDSRL,-EDSCPX) 

EPB»CMPLX(1.0,-0.0) 

GO TO 20 

15 A1=0. 78+0. 03*FREQ-0. 58*1. OE-03* (FREQ* *2) 

A2=0. 97-0. 39*FREQ*1.0E-02+0. 39*1. 0E-03*( FREQ* *2) 

B1=0.31-0.05*FREQ+0.87*1.0E-03*CFREQ*»2) 

A=1.0+(1.83*RHOS) + (0.02*A1*(11V**1.015))+B1 

B=0.073*A1 

C=0.073*A2 

X=1.31 

F0=9.07 

Z=FREQ/F0 

EWSRL=A+C(B*(HV**X))/(1.0+Z*Z)) 
EWSCPX=(C*Z*MV**X)/(1 .0+Z*Z) 
EP2=CMPLX(EySRL,-EUSCPX) 

T=-1.0 

EWINF=4.9 

EW0=88.045-0.4147*T+« .295E-04*T*T+1 .075E-05*T**3 
X1=2.0*PI 

A=1.1109E-10~3.824E-12*T 

B=6 .938E-14*T*T~5 .096E~16*T**3 

T0W=(A+B)/X1 

C-(X*F1(I)*T0W)**2 

Dal.Q+C 

EPURL=EUINF+C(EUO-EUINF)/D) 
EPUCX=(X1*F1(I)*T0U*(EW0-EWINF))/D 
EPU=CMPLX ( EPURL , - EPWCX ) 

A1=(1.0-(3.0*MV/100.0))*(EPW-1.0)-1.0 
EPB=(-Al+CS0RT(Al*Al+8.0*EPU))/4.0 
20 LMDA0=3.0E08/F1(I) 

EPSR«EP2/BPB 

AK=(EPSR-l.d)/(EPSR+2.0) 
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22 

23 


200 

300 

400 

500 


Z»AIMAG(-AK) 

KAI=(6.0*PI*VI*Z*CABS(EPB))/LKDA0 
IFCMV.BQ.O.O)GO TO 22 

KAB=2 . 0*K1*( 1 . 0-VI ) *ABS(AIMAG(CSORT(EPB) ) ) 

GO TO 23 

KAB=0.0 

KA=KAI+KAB 

CA=32.0*PI**4*R**3*VI*(CABS(AK))**2*CABS(EPB)**4 

DA=L1IDA0**4 

KS=CA/DA 

KE=KA+KS 

ALB=KS/KE 

DELTP=1.0/KE 

WRITE(25,*)F1(I),MV,R1,RH0S,KA,KS,KE,DELTP 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

STOP 

END 


C****** 
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APPENDIX - A6 


******************************************************** 
THIS PROGRAM COMPUTES THE BRIGHTNESS TEMPERATURE OF 
TUO-LAYERED MODEL / PROGRAMMED BY SUSHI L K. SRIVASTAV 
******************************************************** 

VARIABLES INVOLVED ARE FREQUENCY, REAL AND IMAGINARY 
PARTS OF DIELECTRIC CONSTANT OF PURE ICE, DENSITY OF 
PURE ICE AND SNOU, RADIUS OF SNOU PARTICLES, LIQUID 
UATER CONTENT OF SNOU, AND TEMPERATURE. 

X************* * * **************************************** 

FI » FREQUENCY IN GHz 

EPICPX = DIELECTRIC LOSS FACTOR OF PURE ICE 

R = RADIUS OF SNOU PARTICLES 

MV = SNOU LIQUID UATER CONTENT 

RHOI = DENSITY OF PURE ICE 

RHOS = DENSITY OF SNOU 

VI = VOLUME FRACTION OF ICE 

EPB = COMPLEX DIELECTRIC CONSTANT OF BACKGROUND MEDIUM 
KAI .= ABSORPTION COEFFICIENT OF ICE PARTICLES 
KAB = ABSORPTION COEFFICIENT OF BACKGROUND MEDIUM 
KS = SCATTERING COEFFICIENT OF THE MEDIUM 
THETO = ANGLE OF INCIDENCE 

EP3 * COMPLEX DIELETRIC CONSTANT OF GROUND MEDIUM 

t******************************************************* 


DIMENSION F1C200),EPICPX(200) 

REAL K1 , K2 , KA , KE , KS , L , MV , LHDAO , KAI , KAB , MVl 
COMPLEX EP2,EP3,R2V,R1V,R2H,R1H,AK,EPSR,EPB,EPU,A1 
OPEN(UNIT=21 ,FILE='SKS6 .DAT’ ) 
OPEN(UNIT=22,FILE=’S18.RES’ ) 
OPEN(UNIT=25,FILE«’S14.RES’ ) 

DO 400 IM=1,2 
DO 300 1=1,35 
READ(21,*)F1(I),EPICPX(I) 

F=F1(I)*1 .0E09 
TI=272.0 
TG=272.0 
Dl=0.5 

D11=D1*100.0 

RHOS=0.25 

UE=DH*RHOS 

RHOI =0.916 

R*l.fl*1.0E“0.1 

PI=3. 142857 

FREQ=F/1 .0E09 

K1=(2.0*PI*F)/3.0E08 

VI=RHOS/RHOI 

EP3=CMPLX(3.0,-0.05) 

IFCMV.GT.O.OJGO TO 15 
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EDSRL* (1.0+0. 508 *RHOS ) * * 3 

BC«EDSRL**2 

CC«2.0*EDSRL 

DC=3.15+CC 

EC=3.15+(2.0*BC) 

EDSCPX«3.0*VI*EPICPXCI)*BC*CCC+1 .0)/CDC*EC) 
EP2=CMPLX(EDSRL. -EDSCPX) 

EPB=CIIPLX(1 .0,-0.0) 

GO TO 20 

15 A1=0. 78+0. 03*FREQ-0. 58*1. 0E-03*(FREQ**2) 

A2=0.97-0.39*FREQ*1.0E-02+0.39*1.0E-03*CFREQ**2) 

B1*0.31-0.05*FREQ+0.87*1.0E-03*(FREQ**2) 

A=1.0+(1.83*RHOS)+(0.02*A1*(MV**1.015))+B1 

B=0.073*A1 

C=0.073*A2 

X=1.31 

F0=--9.07 

Z=FREQ/F0 

EUSRL=A+((B*(MV**X))/(1 .0+Z*Z)) 

EUSCPX= (C*Z*MV**X) / C 1 . 0+Z*Z ) 

EP 2 =CI1PLX ( EUSRL , - EUSCPX ) 

T=-l .0 
EWINF=4.9 

EU0=88.045-0.4147*T+6.295E-04*T*T+1.075E-05*T**3 

X=2.0*PI 

A=1 .1109E- 10-3.82 4E- 12 *T 

B=6. 938E-14*T*T-5.096E-16*T**3 

TOU=(A+B)/X 

C=CX*F*T0W)**2 

D=1.0+C 

EP«RL=E«INF+C(EUO-EWINF)/D) 

EPWCX=(X*F*TOW*(EWO-EWINF))/D 

EPW=CI!PLX(EPWRL,-EPWCX) 

A1=(1.0-(3.0*1IV/100 .0))*(EPW-1.0)~1.0 
EPB=(-Al+CSQRT(Al*Al+8.0*EPU))/4.0 
20 LMDA0=3.0E08/F 

EPSR=EP2/EPB 
AK=(EPSR-1.0)/CEPSR+2 .0) 

Z=AIMAG(-AK) 

KAI=,(6.0*PI*VI*Z*CABSCEPB))/L11DA0 
IF(MV.EQ.O.O)GO TO 22 

KAB=2 . 0*K1 * ( 1 . O-VI ) *ABS(AIMAGCCSQRT(EPB) ) ) 

GO TO 23 

22 KAB=0.0 

23 KA=KAI+KAB 

CA=32.0*PI**4*R**3*VI*(CABS(AK))**2*CABSCEPB)**4 

DA=LMDA0**4 

KS=CA/DA 
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KE«KA+KS 

ALB=KS/KE 

THET0®0.0 

30 THET1=THET0*PI/180.0 

THET2=ASINCCABS(CS0RT(1.0/EP2))*SIN(THET1)) 

THET3»ASIN(GABSCCSQRTCEP2/EP3))*SIN(THET2)) 

L=EXP(KE*D1*1 . 0/COS(THET2 ) ) 

AmiO=4.O*PI*l.OE-07 

EP0*8.854E-12 

Z1=SQRT(AMU0/EP0) 

Z1V=Z1*C0S(THET1) 

Z2V=(Z1*C0S(THET2)/CSQRT(EP2)) 

Z3V= (Z1 *C0S(THET3 ) /CSQRT(EP3 ) ) 

R2V=(Z2V-Z3V)/(Z2V+Z3V) 

R1V=(Z1V-Z2V)/(Z1V+Z2V) 

REFL2V=(CABSCR2V))**2 

REFL1V=CCABS(R1V))**2 

AA1=1 -O-REFLIV 

BBl=i.Q 

CCl=i .0+CREFL2V/L) 

DD1=1 .0-(l .0/L) 

EB1=C1,0-ALB) 

FF1=((1.0-REFL2V)/L) 

EMISV=AA1*((CC1*DD1*EE1*TI)4FF1*TG)/BB1 

Z1H=Z1/C0SCTHET1) 

Z2H«(Z1/(C0SCTHET2)*CSQRT(EP2))) 

Z3H=(Z1/(C0SCTHET3)*CSQRT(BP3))) 

R2H=(Z3H-Z2H)/(Z3H+Z2H) 

R1H=CZ2H-Z1H)/(Z2H+Z1H) 

REFL2H=(CABSCR2H))**2 

REFL1H=CCABSCR1H))**2 

AA2=1 .0-REFLlH 

BB2:=1.0 

CC2=1 .0+(REFL2H/L) 

DD2=1.0-(1.0/L) 

EE2=(1.0-ALB) 

FF2=C(1.0-REFL2H)/L) 

EMISH=AA2*(CCC2*DD2*EE2*TI)4FF2*TG)/BB2 
TBVMTBH = EM I S V- EM I SH 

WRITE(22,*)IM,F1(I),THET0,EMISV,ENISH 
THETb=THET0+4.0 
if(theto.le.&o.o)go to 30 
300 CONTINUE 

400 CONTINUE 

STOP 
END 


Q* * j(r * * * 
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APPENDIX - A7 


c***** this program computes the brightness temperature of 
c***** multi-layered models/ programmed by sushi L K. SRIVASTAV- 

Q* * t It t * 1 1 1 * * l! K * * * * * It It * * 1c t ic It t * it * * Hie im ic * it * * * Hi * * * ** it * * * * * * * it * * * * t * * 

c***** VARIABLES ARE SAME AS GIVEN IN THE PROGRAM FOR TWO- 
C***** LAYERED MODEL. SIG IS THE SURFACE HEIGHT VARIATION. 

Qfet 1 1 1 1 1 1 1 * 1 1 1 * t it t 1 1 1 1 1 1 1 1 1 ic -k ic t * 1 1 ie 1 1 1 ict t ttic * * t * * t * icictic t icic 

DIMENSION F1(150),F(150),EPICPX(150),RHOS(150,10) 

DIMENSION R(150,10),H(150,10),EP(150,10),Z(i50,10) 

DIMENSION KE(150,10),ALB(150,10),THET(150,10),X(150,10) 
DIMENSION REFLV(150,10),ENISV(150,10),RR1(150,10) 

DIMENSION EMIS1(150),TB1(150),KA(150,10),KS(150,10) 

DIMENSION L(150,10),ZV(150,10),RFL(10),Ri(150,10) 

DIMENSION RR2(150,10),ZH(150,10),REFLH(150,10),y(150,10) 
DIMENSION EMISH(150,10),EMIS2(150),TB2C150),MV(150,10) 
DIMENSION SIG(150,10),REF1H(150,10),RBF1V(150,10),KAB(150,10) 
DIMENSION K(150,10),HDASH(150,10),FACT(150,10),KAI(150,10) 
DIMENSION MV1(150,10) ,EPW(150,10),BPB(150,10),EPSRC150,10) 
REAL L,K,KA,KS,KE,KAI,LMDAO,TB,MV,MVl,KAB 
COMPLEX EP , AK , ZV , ZH , AB , EPB , EPU , EPSR 
OPEN(UNIT=21,FILE=’SKS6.DAT') 

OPEN(UNIT=22, FILE* •SKS7.DAT’) 

OPENCUNIT=24,FILE*’S26.RES' ) 

DO 700 1=1,115 
READ(21, *)F1(I) ,EPICPX(I) 

DO 600 J=2,5 

READ(22 ,*)RHOSCI , J) ,R1(I , , J) 

600 CONTINUE 

700 CONTINUE 

PI=3. 142857 

EP0=8.854E-12 

AMU0=4.0*PI*1.0E-07 

T=272.0 

C=3.0E08 

RHOI=0. 916 

DO 550 1=1,1 

EP(I ,l)=CMPLX(i.0,-0.0) 

EP(I ,6)=CMPLX(3.5,-0.05) 

Z( I , 1 )=S0RT(AMU0/EP0) 

550 CONTINUE 

DO 545 1=1,115 
DO 540 J=2,5 
SIGCI ,1)=0.01 
SIG(I,J)=0.0 
540 CONTINUE 

545 CONTINUE 

DO 481 1=1,1 
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DO 480 J=2,5 
I!V(I,J) = 1.0 

480 CONTINUE 

481 CONTINUE 

DO 500 1=1,1 

DO 400 J=2,5 

VI=RKOSCI , J)/RHOI 

FCI)=F1(I)*1,0E09 

K(I,1)=2.0*PI*F(I)/C 

IF(MV’(I,J).NE.O.O)GO TO 450 

BDSEL=(1.0+0.508*RHOS(I , 

BC=EDSRL*EDSRL 

CC=2.0*EDSRL 

DC=3.15+CC 

EC=3.15+(2.0*BC) 

EDSCPX=3.0*VI*BPICPX(I)*BC*(CC+1. )/(DC*EC) 

EP(I , J)=CMPLX(EDSRL,-EDSCPX) 

GO TO 400 

450 AI1=0.78+0.03*F1(I)-0.58*1.0E-03*CF1(I)**2) 

AI2=0. 97-0.39*Fl(I)*l .0E-02+0.39*1.0E-03*CFl(I)**2) 

BI1=0. 31-0,05*F1(I)+0.87*1.0E-03*(F1(I)**2) 

AI3=1.0+C1.83*RHOS(I, J))+(0.02*AIi*(MV(I , J)**1.0i5))+BIl 

BI2=0.073*AI1 

XI1=0.073*AI2 

XI=1.31 

F0I=9.07 

AZI=F1(I)/F0I 

EDSRL=AI3+((B12*(MV(I ,J)**XI))/(1 .0+AZI*AZI)) 
EDSCPX=CXI1*AZI*I!V(I , J)**XI)/(1.0+AZI*AZI) 

EP(I , J)=CWPLX(EDSRL,-BDSCPX) 

URITE(24,*)E?(I,J) 

00 CONTINUE 

DO 396 J=2,6 

K(l , J)=KCI ,1)*CABS(CSQRT(EP(1 , J))) 

396 CONTINUE 

DO 390 J=2,5 

IF(MVCI,J).EQ.O.O)GO TO 395 

Tl=T-273.0 

EUINF=4.9 

BU0=88.045-0.4147*T1+6.295E-04*T1*T1+1.075E-05*T1**3 

XU=2.0*PI 

AU=1 . 1109E-i0-3 . 824E-12*T1 

BU=6. 938E“14*Tl*Tl-5.096E~16*Ti**3 

TOW=(AW+BW)/Xy 

CW=(XU*F(I)*TOW)**2 

DW=1.0+CW 

EPyRL=EWINF+(CEWO-£MINF)/DU) 

EPWCPX=CXU*FCI'5*TOW*(E«0-EUINF))/DW 
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EPWCI, J)=CMPLX(EPWRL,-EPWCPX) 

MVi(:: , j)=MV(i , j)/ioo.o 

AB=(1.0-3.0*MV1(I, J))*(EPU(I, J)--1.0)-1.0 
EPB(I ,J)=(-AB+CS0RT(AB*AB+8.0*EPU(I , J)))/4.0 
GO TO 390 

395 EPB(I , J)=CMPLX(1.0,-0.0) 

390 CONTINUE 

DO 300 J=2,5 
VI=RHOS(I, J)/RHOI 
R(I, J)=R1(I, J)*1.0E-03 
LMDAO=C/F(I) 

EPSRCI , J)=EPCI , J)/BPB(I, J) 

AK=(BPSR(I , J)-1.0)/(EPSR(I, J)+2.0) 

Z2=AIMA6(-AK) 

KAI(I , J)=(« .0*PI*VI*ZZ*CABS(EPB(I , J)))/LMDAO 
IFCMVCI, J).EQ.O.O)GO TO 380 

KABCI , J)=2.0*K(I, J)*(1.0-VI)*ABS(AIMAG(CSQRT(EPB(I, J)))) 

GO TO 370 

380 KABCI,J)»0.0 

370 KA(I , J)=KAI(I , J)+KAB(I , J) 

CA=32 .0*PI**4*R(I , J)**3*VI*CCABS(AK))*’'2*CABS(EPB(I , 

DA=LHDA0**4 

KS(i , J)=CA/DA 

KECI , J)=KS(I, J)+KA(I , J) 

ALBCi; , J)=KSCI , J)/KE(I , J) 

300 CONTINUE 

THET0=0.0 

250 THETCI ,1)=THET0*PI/180.0 

DO 200 J=2,6,l 

X(I , J)=CABSCCSORT(EP(I , ( J-1 ) ) /EP( I , J) ) ) 

Y(I , J)=XCI , J)*SIN(THET(I , (J-1))) 

THET(I , J)=ASIN(Y(I , J)) 

200 CONTINUE 

DO 150 J=2,5 

L(I , J)=EXP(KECI , J)*H(I , J)/COSCTHET(I , J))) 

C WRITE(24,*)LCI , J) 

150 CONTINUE 

ZVCI ,1)=Z(I,1)/C0S(THET(I .1)) 

ZHC I , 1)=ZC I , 1)*C0S(THETC I , 1) ) 

DO 75 J=2,6 

ZV(I , J)=(ZV(I i (J-1))*C0S(THET(I , J))/CSQRT(EP(I , J))) 

2HCI , J)=CZH(I , (J-1))/(C0S(THETCI , J))*CSQRTCEPCI, J)))) 

C WRITE(24,*)ZV(I,J),ZH(I ,J) 

75 CONTINUE 

DO 50 J=l,5 

REFIVCI . J)=CABSC(2V(I. J)-2V( I, (J+1)))/CZV(I , J)+ZV(I , CJ + 1))))**2 
REFIHCI , J)=CABSCCZH(I, (J + 1))-ZH(I , J))/(ZH(I , J)+ZH(I, (J + l))))**2 
50 CONTINUE 

DO 35 J=l,5 
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KDASHCI, J) = 4.0*K(I , J)*KCI, J)*SIG(I , J)*SIGCI, J) 
FACTCI , J)-EXPC-HDASHCI , J)*COSCTHET(I , 

REFLVCI , J)=REF1V(I , J)*FACT(I, J) 

REFLHCI , J)=REF1H(I , J)*FACT(I , J) 

35 CONTINUE 

DO 40 M=l,4 
J=6-M 

AAl=i . 0-REFLVC I , ( J-1 ) ) 

AA2“1.0-EEFLH(I, (J-l)) 

BBls^i^O-CREFLVd , J)*REFLV(I, (J~1))/(L(I , J)*LCI , J))) 
BB2=1.0-CREFLH(I , J)*REFLH(I,(J-1))/(L(I, J)*LCI, J))) 
CC1=1.0+<REFLV(I , J)/LCI, J)) 

CC2 = i.0 + (REFLH(I, J)/L(Id)) 

DDl"1.0-(l-0/L(I.J)) 

EE1=1.0~ALB(I, J) 

FF1=((1.0-REFLV(I, J))/L(I , J)) 

FF2=((i.0-REFLH(I, J))/LCI, J)) 

EMISVCI , J)=AA1*((CC1*DD1*EE1)+FF1)/BB1 
EMISHCI , J)=AA2*CCCC2*DD1*EE1)+FF2)/BB2 
REFLV(I ,CJ-1))=1.0-EMISV(I,J) 

REFLHCI , (J-1))=1.0-BMISH(I , J) 

40 CONTINUE 

EMISl(I)=1.0-REFLHCI,i) 

EMIS2 ( I )=1 . 0-REFLV( 1,1) 

TB1(I)=EMIS1(I)*T 

TB2CI)=EI!IS2CI)*T 

URITE(24,*)F1<I) ,THET0,TB1(I) ,TB2(I) 
IFCTHET0,GT.50.0)GO TO 500 
THET0=THET0+4.0 
GO TO 250 
500 CONTINUE 

STOP 
END 



